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1 Introduction 

In this paper we use the method of symmetry reduction for exterior differ- 
ential systems to obtain a far-reaching generalization of Vessiot's integration 
method [57], [SB] for Darboux integrable, partial differential equations. This 
group-theoretic approach provides deep insights into this classical method; un- 
covers the fundamental geometric invariants of Darboux integrable systems; 
provides for their algorithmic integration; and has applications well beyond 
those currently found in the literature. In particular, our integration method 
is applicable to systems of hyperbolic PDE such as the Toda lattice equations, 
[T9] , [20] , [22] , 2-dimensional wave maps [3] and systems of overdetermined PDE 
such as those studied by Cartan [5J. 

Central to our generalization of Vessiot's work is the novel concept of a 
superposition formula for an exterior differential system X on a manifold M. 
A superposition formula for I is a pair of differential systems W, W, defined on 
manifolds M\ and M2, and a mapping 

E:MixM 2 ->M (1.1) 

such that 

£*(X) C 7T*(W) +7r|(VV). (1.2) 

Here irl (W + (W) is the differential system generated by the pullbacks of W 
and >V to the product manifold Mi x Mi by the canonical projection maps tt\ 
and 7T2. It is then clear that if cf>: N\ — > Mi and cf>: N2 — > M2 are integral 
manifolds for W and W, then 

= So $)-. N r X iV 2 -> M (1.3) 

is an (possibly non-immersed) integral manifold of X. 
In this article we shall 
[i] establish general sufficiency conditions (in terms of geometric invariants of 
the differential system I) for the existence of a superposition formula; 

[ii] establish general sufficiency conditions under which the superposition for- 
mula gives all local integral manifolds of X in terms of the integral manifolds 
of 7T* (W) +7T2 (W) on Mi x M 2 (in which case we say that the superposition 
formula is surjective); 

[iii] provide an algorithmic procedure for finding the superposition formula; and 
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[iv] demonstrate the effectiveness of our approach with an extensive number of 
examples and applications. 
Differential systems admitting superposition formula are easily constructed 
by symmetry reduction. To briefly describe this construction, let G be a sym- 
metry group of a differential system W on a manifold N. We assume that the 
quotient space M = N/ G of N by the orbits of G has a smooth manifold struc- 
ture for which the projection map q: N — > M is smooth. We then define the 
G reduction of W or quotient of W by G as the differential system on M 
given by 

W/G = {weff(M)|q»6W}. (1.4) 

The traditional application of symmetry reduction has been to integrate W by 
integrating W/G. See, for example, [T]. 

But now suppose that differential systems W and W on manifolds M\ and 
M2 have a common symmetry group G. Define the differential system VV = 
7r* (W) + vrj (W) on Mi x M 2 and let G act on Ml x M 2 by the diagonal action. 
Then, by definition, the quotient map q: Mi x M 2 — > M = (Mi x M 2 )/G 
defines a superposition formula for the quotient differential system W/G on M. 
In this paper we discover the means by which the inverse process to symmetry 
reduction is possible, that is, we show how certain general classes of differential 
systems X on M can be identified with a quotient system W/G in which case 
the integral manifolds of I can then be found from those ofW and W . 

To intrinsically describe the class of differential systems for which we shall 
construct superposition formulas, we first introduce the definition of a decom- 
posable differential system. 

Definition 1.1. An exterior differential system X on M is decomposable of 

type [p, q], where p, q > 2, if about each point x € M there is a coframe 

e\ e r , o\ a\ 59, (1.5) 

such that X is algebraically generated by 1-forms and 2-forms 

x = {~e\ e r , ft 1 , ti s , h\ ft}, (1.6) 

where s,t>l,£l a £ fl 2 (a\ . . . ,a p ), and fl a G fi 2 ^ 1 , . . .,a q ). The differential 
systems algebraically generated by 

V = {6\ a a ,n a } and V = {9\a a ,fl a } (1.7) 
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are called the associated singular differential systems for 2 with respect to 
the decomposition (|1.6p . 

With the goal of constructing superposition formulas, we have found it most 
natural to focus on the case where 1 is decomposable (but not necessarily Pfaf- 
fian) and V and V are (constant rank) PfaffiarQ All of the examples we consider 
are of this type. Note that any class r hyperbolic differential system, as defined 
in [6], is a decomposable differential system and that the associated character- 
istic Pfaffian systems coincide, for r > 0, with the singular systems ( |1.7[) . 

The definition of a Darboux integrable, decomposable differential system is 
given in terms of its singular systems. For any Pfaffian system V, let V^ 00 ' 
denoted the largest integrable subbundle of V. The rank of V^ 00 ^ gives the 
number of functionally independent first integrals for V. By definition, a scalar 
second order partial differential equation in the plane is Darboux integrable if 

A v 

the associated singular Pfaffian systems V and V each admit at least 2 (func- 
tionally independent) first integrals. Thus, in order to generalize the definition 
of Darboux integrablity, we must determine the required number of function- 
ally independent first integrals necessary to integrate a general decomposable 
Pfaffian system. We do this with the following definitions. 

Definition 1.2. A pair of Pfaffian systems V and V define a Darboux pair 

if the following conditions hold. 

[i] V + V {co) =T*M and V + V (oo) =T*M. (1.8) 

[ii] f W n^ = {0}. (1.9) 

[iii] duj e n 2 (V) + fl 2 {V) for all lo efl^VnV) . (1.10) 

Definition 1.3. Let 2 be a decomposable differential system and assume that 

/\ — 

the associated singular systems V and V are Pfaffian. Then X is said to be 
Darboux integrable if {V, V} define a Darboux pair. 

Property [i] of Definition 11.21 is the critical one - it will insure that there 
are a sufficient number of first integrals to construct a superposition formula. 
Property [ii] is a technical condition which states simply that V and V share no 

1 We use the term Pfaffian system to designate either a constant rank subbundle V of T* M 
or the differential system, denoted by the corresponding calligraphic letter V, generated by 
the sections of V. 
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common integrals - this condition can always be satisfied by restricting V and 
V to a level set of any common integrals. The form of the structure equations 
for V n V required by property [iii] is always satisfied when V and V are the 
singular Pfafhan systems for a decomposable differential system T. 
Our main result can now be stated. 

Theorem 1.4. LetX be a decomposable differential system on M whose asso- 
dated singular Pfaffian systems { V, V } define a Darboux pair. Then there are 
Pfaffian systems W\ and W2 on manifolds M 1 and M 2 which admit a common 
Lie group G of symmetries and such that 

[i] the manifold M can be identified (at least locally) as the quotient of M\ x M2 
by the diagonal action of the group G; 

[ii] 1= (7r 1 (W + 7r*(VV))/G; and (1.11) 

[iii] the quotient map q: Mi x M2 — > M defines a surjective superposition for- 
mula for X. 

The manifolds M\ and M2 in Theorem 1 1.41 are simply any maximal integral 
manifolds for V^°°^ and y^ 00 ^ and the Pfaffian systems W\ and W2 are just the 
restrictions of V and V to these manifolds. But the proper identification of the 
Lie group G and its action on M x and M 2 is not so easy to uncover. This is done 
through a sequence of non-trivial coframe adaptations and represents the prin- 
ciple technical achievement of the paper (See Theorem 14. II and Definition 15. 7[) 
We call G the Vessiot group for the Darboux integrable, differential system 
1. 

The paper is organized as follows. In Section 2 we obtain some simple 
sufficiency conditions for a differential system to be decomposable and we give 
necessary and sufficient conditions for a Pfaffian system to be decomposable 
(Theorem 12.31) . We also introduce the initial adapted coframes for a Darboux 
pair (Theorem l2.9p . In Section 3 we answer the question of when the symmetry 
reduction of a Darboux pair is also a Darboux pair (Theorem 13. 2p and we 
use this result to give a general method for constructing Darboux integrable 
differential systems (Theorem l3.3|) . Section 4 establishes the sequence of coframe 
adaptations leading to the definition of the Vessiot group G and its action on 
M and hence on Mi and M2 (Theorem 14. ip . In Section 5 we construct the 
superposition formula (Thcorcm l5.10p and prove that it may be identified with 
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the quotient for the diagonal action of the Vessiot group. This proves Theorem 
PI 

A substantial portion of the paper is devoted to examples. These exam- 
ples, given in Section 6, provide closed-form general solutions to a wide range 
of differential equations, explicitly illustrate the various coframe adaptations 
of Section 4, and underscore the importance of the Vessiot group as the fun- 
damental invariant for Darboux integrable differential systems. The first two 
examples are taken from the classical literature. The equation considered in 
Example 16. II is specifically chosen to illustrate in complete detail all the various 
coframe adaptations used to determine the superposition formula. In Example 
I6.2l we integrate all non- linear Darboux equations of the type u xx = f{u yy ) (with 
second order invariants) and relate the integration of these equations to Car- 
tan's classification of rank 3 Pfaffian systems in 5 variables [7] . This connection 
between the method of Darboux and Cartan's classification is apparent from 
our interpretation of the method of Darboux in terms of symmetry reduction of 
differential systems. 

In Examples 16.31 and 16.41 we present a number of examples of Darboux in- 
tegrable equations where the unknown function takes values in a group or in a 
non-commutative algebra. Example [63] provides us with a system whose Vessiot 
group is an arbitrary Lie group G. 

Some of the simplest examples of systems of Darboux integrable partial dif- 
ferential equations can be constructed by the coupling of a nonlinear Darboux 
integrable scalar equation to a linear or Moutard-type equation. These are pre- 
sented in Example 16.51 It is noteworthy that for these equations the Vessiot 
group is a semi-direct product of the Vessiot group for the non-linear equa- 
tion with an Abelian group. The representation theoretical implications of this 
observation will be further explored elsewhere. 

In Example 16.61 we illustrate the computational power of our methods by 
explicitly integrating the B2 Toda lattice system. While the general solution 
to the A n Toda systems were known to Darboux, this is, to the best of our 
knowledge, the first time explicit general solutions to other Toda systems have 
been given. Based on this work we conjecture that the Vessiot group for the q 
Toda lattice system, where g is any semi-simple Lie group, is the associated Lie 
group G itself. In Example 16.71 we integrate a wave map system - this system 
admits a number of interesting geometric properties which will be explored 
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in detail in a subsequent paper. Finally, in Example 16.81 we present some 
original examples of Darboux integrable, non-linear, over-determined systems 
in 3 independent variables. 

The extensive computations required by all these examples were done using 
the DifferentialGeometry package in Maple 11. This research was supported by 
NSF grant DMS-0410373 and DMS-0713830. 
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2 Preliminaries 

2.1 Decomposable Differential Systems 

In this section we give some simple necessary conditions, in terms of the 
notion of singular vectors, for a differential system^) 2" to be decomposable (see 
Definition II. 1| ; we give sufficient conditions for a Pfaffian system to be decom- 
posable; and we address the problem (see Theorem 12. 6|) of re-constructing a 
decomposable differential system from its associated singular systems. 

Let X be a differential system on M. Fix a point x in M and let 

El(l) = { X e T X M | 6(X) = for all 1-forms 6 el}. 

The polar equations determined by a non-zero vector X e E*(l) are, by 
definition, the linear system of equations for Y e E%(1) given by 

0{X, Y) = for all 2-forms 9 el. 

Then X e E^(l) is said to be regular if the rank of its polar equations is 
maximal and singular otherwise. 

Lemma 2.1. Let 1 be a decomposable differential system of type [p,q\. Then 
E\(l) decomposes into a direct sum of p and q dimensional subspaces 

El(l)=S 1 ®S 2 (2.1) 

such that 

[i] every vector in Si and every vector in S 2 is singular; and 

[ii] every 2 plane spanned by any pair of vectors X e Si and Y e S 2 is an 

integral 2-plane for 1. 

Proof. Let { , , da a } denote the dual frame to (II. 5p . Then E\ (1) = 
span{ d^i , da <* } and the polar equations for X — X\ + X 2 , where Xi — t a d a ° 
and X 2 = s a d a o, , are 

Q a (Xi,Y)^0 and n a {X 2l Y) = 0. 

Then, clearly, the vectors X\ and X 2 are singular vectors and, moreover, the 
plane spanned by X\ 1 X 2 is an integral 2-plane. | 

2 We assume that all exterior differential systems are of constant rank. 
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Remark 2.2. The differential system V denned by (|1.7[) is the smallest differ- 
ential system containing X and Ann(S , 2) and it is for this reason that we have 
opted to called V (and V) the singular differential systems associated to the 
decomposition of X. 

The necessary conditions for decomposability given in Lemma [2.1l are seldom 
sufficient. However, if I is a Pfaffian system and properties [i] and [ii] of Lemma 
12. II are satisfied, then there exists a local coframe 

9\...,6 r , a\...,a^a\...^ 

on M with / = { 6 l } and with structure equations 

dQ i = A\ b a a Aa b + B^ a a A o (i mod /. (2.2) 

At each point xq of M define linear maps 

A=[Ai b (x )}:R r ^A 2 (Rf) and B = [B' l af3 (x Q )} : R r -> A 2 (R 9 ). (2.3) 

Theorem 2.3. Let I be a Pfaffian system. Then I is decomposable if properties 
[i] and [ii] of Lemma \2.1\ are satisfied; the matrices A l ab and B l a a given by (12. 2|) . 
are non-zero, constant rank; and 

dim ( ker(A) + ker(B)) = rank /. (2.4) 

Proof. Equation (|2.4|) implies that we can find r linearly independent r dimen- 
sional column vectors T 1 ,. . . , T r \. . . , T r2 , . . . , T r such that 

kei(A) n ker(S) = span{ T 1 , ... , T ri }, 
ker(S) = span{ T 1 , . . . , T r \T ri+1 , ...,T r2 }, and 
ker(A) = span{ T 1 , . . . , T ri , T r2+1 , . . . , T r }. 

The 1-forms 9 l = t) 9 j , where T l = [tj], then satisfy 

d9 l = forz = l...n, 

d¥ = t'jAl b a a Aa b for i = n + 1 . . . r 2 , (2.5) 
=t)B J a0 a a Act' 3 for i = r 2 + 1 . . . r, 

and the decomposability of / follows by taking fl l — tjA J ab a a A a b for i — 
r\ + 1 . . . r 2 and (l l — t l jB 3 a/3 a a A a h for i = r 2 + 1 . . . r. That the integers s and 
t in Definition 11.11 satisfy s,t > 1 follows from the fact that and B l a ^ are 
non-zero and that [tj] is invertible. | 
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Corollary 2.4. IfX is a decomposable Pfaffian system, then the singular sys- 

A '■■ 

terns V and V are also Pfaffian. 

Proof. This is immediate from (|1.7[) and (|2.5[) . | 

To each decomposable exterior differential system X we have associated a 

pair of differential systems V and V and, when these are Pfaffian, we have 

a ^ 

defined what it means for { V, V } to define a Darboux pair. In our subsequent 
analysis, the Darboux pair { V, V } will be taken as the fundamental object of 
study. From this viewpoint, it becomes important to address the problem of 
reconstructing the EDS X from its singular systems. To this end, we introduce 
the following novel construction. 

Definition 2.5. Let V and V be two differential systems on a manifold M. 

Define a new differential system K, = V ffl V to be the exterior differential system 

a 

whose integral elements are precisely the integral elements of V , the integral 
elements of V , and the sum of integral elements of V with integral elements of 
V. 

This definition is motivated by two observations. First we note that if X is a 

a 

decomposable exterior differential system with singular Pfaffian systems V and 

V then in general X properly contains the Pfaffian system with generators VD V 
and is properly contained in the EDS Vn V so that neither of these set-theoretic 
constructions reproduce X. 

Secondly, one can interprete the EDS K, in Definition 12.51 as satisfying an 
infinitesimal superposition principle with respect to { V, V } at the level 
of integral elements. Of course, this by itself does not imply that K. admits a 
superposition principle for its integral manifolds. The main results of this article 
can then be rc-formulated as follows: if the EDS K. admits a superposition 
principle for its integral elements with respect to a Darboux pair { V, V}, then 
K, admits a superposition principle for its integral manifolds. 

Theorem 2.6. IfX is a decomposable exterior differential system with singular 

A ^ A ^ 

differential systems V and V, then X = V Ffl V. 

Proof. We first remark that the differential system V ffl V is always contained in 

V ffl V. Let {d 1 , a a , a a } be the coframe given by Definition ll.il By definition, 

A ^ 

the algebraic generators for V and V are 

V = {0\ a a , fV 3 } and V = {§, a a , &} 
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from which it is not difficult to show that 

A 

Since <9jc* defines a 1-integral clement for V and da<* defines a 1-integral element 
for V, the 2-plane spanned by {c%o, 9^°} must be a integral element for V ffl V 
and therefore the forms a a A a a ^ V FR V. In view of (|1.6[) . we conclude that 

To complete the proof, it suffices to check that every (r + s)-form u E X 
vanishes on every plane E r + E s , where E r is an r- integral plane of V and E s 
is an s-integral plane of V. The form lj is a linear combination of the wedge 
product of the 1-forms in (|1.6[) with arbitrary (7' + s — l)-forms p and the wedge 
product of the 2-forms in (|1.6[) with arbitrary (r + s — 2)-forms r. It is clear 

~ , a. , 

that (6?* A p)\(E r +E 3 ) ~ 0- The va hies of (ST A ^k^+b^ can be calculated as 
linear combinations of terms involving 

A. 1 A. A. At A* "V , -v" /S, /\ , v v 

r! b (x,r), fi b (x,r), fi 6 (x,n ft b (x,n 

where X, F 6 E r and X, F G E s . These terms all vanish and in this way we 
may complete the proof of the theorem. | 

2.2 The first adapted coframes for a Darboux pair 

Let { V, V } be a Darboux pair on a manifold M. In addition to the prop- 
erties listed in Definition \1.2\ we shall assume that the derived systems 
and V { °°\ as well as the intersections V [oc) H V, V H V r( °° ) and V" H 7 are all 
(constant rank) subbundles of T*M. 

To define our initial adapted coframe for the Darboux pair { V, V }, we first 
chose tuples of independent 1-forms fj and f] which satisfy 

V (oo) n V = span{f)} and t> H V (co) = spanjj)}. (2.6a) 

Property (| 1 . 9|) implies that the forms {77, fj } are linearly independent. Then 
chose tuples of 1-forms 0, <x and <r such that 

V {oo) = span{ <r, 17 }, V (oo) = span{ <r, 17 }, t> n V" = span{ 6, r), rj } (2.6b) 

and such that the sets of 1-forms { <r, f) }, { <r, 17 } and {0, 17, 17} are linearly 
independent. 
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Lemma 2.7. Let { V, V } be a Darboux pair on a manifold M . Then the 1- 
forms {9, & , fj, a, 77} satisfying (|2.6[) define a local coframe for T*M and 

V = span{ 9, <x, fj, fj } and V — span{ 9, 17, a, fj }. (2.6c) 

Proof. From the definitions (|2.6ap and (|2.6b[) it is clear that span{ 9, a, fj, fj } C 

A v 

V and span{ 9, fj, a, f)} C V. Property (|1.8[) then implies that the forms 
{9, <x, 17, <t, 17} span T*M. 

To prove the first of (|2.6c[) , let u; — a.9 + b & + cfj + d& + efj. If cj € V, 
then the 1-forms da eV and therefore, by ((23b)) . d<r € VnV (oo) . Since the 
forms fj and <x are independent, we must have d = and (|2.6c|) is established. 

To prove the independence of the forms { 9, a, fj, rr, fj }, set uj = 0. Then 

A. v 

the argument just given implies that d = so that b<r £ V (IV. Therefore 
b<r 6 V^°°' n V and this forces b = 0. The independence of the forms {9, fj, fj} 
then gives a = c = e = and the lemma is proved. | 

Any local coframe { 9, cr, fj, a, fj} satisfying (|2.6|) is called a O-adapted 

A v 

coframe for the Darboux pair {V, V}. 

We shall also need the definition of a 0-adapted frame as the dual of the 
0-adapted coframe. Let 

TT T> TT TV Tt(°°) T>(°°) Tt(°°) Tr(oo) 

if = ann V, if = ann K , if = ann K , iT = ann K , 

(2.7) 

and K = H (oc) (1 H {oo) . 

The definition of Darboux integrability can be re-formulated in terms of the 
distributions H,H ([H]). If we introduce the dual basis { d e , d~, d~, d^} 
to the 0-adapted coframe {9, ar, fj, a, fj}, that is, 

9(d e ) = l, = 1, V7(^) = l, = 1, W) = l, (2-8) 

with all others pairings yielding 0, then it readily follows that 

H = span{9^}, H (oo) = span{ d g , , }, 
H = span{d-}, ff (oo) = span{ 8 g , 8 & , 9- }, (2.9) 
and K = span{ d e } . 

Any local frame {d e , d~, d~, d^, d^} satisfying (|2.9[) is called a 0-adapted 
frame. 

We shall use the following lemma repeatedly. 
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Lemma 2.8. Let f be a real-valued function on M. If X(f) — for all vector 
fields in H, then df £ V . Likewise, if X(f) = for all vector fields in H, 
then df G V {co) . 

Proof. It suffices to note that if X(f) = for all vector fields X G H, then 
Z(f) = Z-> df = for all vector fields Z G H (oc) and therefore df G ann if (oo) = 
f (oo) . | 

Real- valued functions on M with df G V (or are called first integrals 

for V. We denote the algebra of all first integrals for V by Int(V). 

Our 1-adapted coframe for a Darboux pair is easily constructed from com- 
plete sets of functionally independent first integrals for y*- 00 - 1 and y(°°' > i that is, 
locally defined functions { I } and { I } such that 

V {oo) = span{ d J } and V (co) = span{ dl }. 

Complete the 1-forms { fj } in (|2.6aj) to the local basis (|2.6bp for j/^ 00 ' using 
forms 

<r = dJi (2.10) 

chosen from the { dl}. Next let { 1 2 } be a complementary set of invariants to 
the set { Ji }, that is, { 1} = { t±, I2 }■ Because the forms 77 belong to y°° we 
can write 

r) = Radix + S di 2 . 

Since the 1-forms 17 are independent of the 1-forms 6", the coefficient matrix So 
must be invertible and we can therefore adjust our local basis of 1-forms {fj} 
for V (oo) n V by setting 

fj = di 2 +RdJi = di 2 +R<x. (2.11) 

The exterior derivatives of these forms are 

dfj = dR A a = d^(R) a A a -\ 

and therefore, on account of (jl.lOp , the fact that fj G V^ 00 ^ n V", and (|2.6cp . we 
must have ^(R) = 0. Lemma [2T81 implies that d(R) € anc i consequently 

dfj = Efj A<r + F& Act. (2.12) 

Similar arguments are used to modify the forms ct and fj. These adaptations 
are summarized in the following theorem. 
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Theorem 2.9. Let { V, V} be a Darboux pair on a manifold M. Then about 
each point of M there exists a 0-adapted coframe { 6, <x, fj, <r, fj } with structure 
equations 

da = 0, dfj = E fj A & + Fa A a, 

( 2 - 13 ) 

da = 0, dfj = E fj A a + F a A a, 
d9 = Aa Aa + Ba Aa mod {6, fj, fj}. (2.14) 

We remark that the structure equations for the 1-forms £ VC\V are a conse- 
quence of (jl.lOp and the properties (|2.6p of a 0-adapted coframe. A 1-adapted 
coframe for the Darboux pair { V , V } is a 0-adapted coframe { 9, & : fj, a, fj} 
satisfying the structure equations (12. 13[) and (|2.14l) . 

Corollary 2.10. If{V, V} define a Darboux pair, then V ffl V is always a de- 
composable exterior differential system for which V and V are singular Pfaffian 
systems. 

Proof. It follows from Theorem 12.91 that the algebraic generators for V and V 
are 

A - 

V = {8, fj, fj, a-, Fa A a, B&Acr} and 

V = {6, fj, fj, rr, F&A&, AaAa}. 

From these equations it is not difficult to argue that the differential system VnV 
has generators 

VnV = {0, fj, fj, &A&, F&Afr, Act Act, F&A&, B&A&). 

A repetition of the arguments used in the proof of Theorem 12.61 shows that 

VffiV = {9, fj, fj, F&A&, A&A&, F&A&, B&A&). (2.15) 

It is then clear that V FR V is a decomposable differential system (Definition [T7T]) 
and that V and V are singular Pfaffian systems. | 

Example 2.11. There are many examples of differential equations which can 
be described either by Pfaffian differential systems or by differential systems 
generated by 1-forms and 2-forms. Our definition of Darboux integrability in 
terms of decomposable differential system is such that theseequations will be 
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Darboux integrable irrespective of their formulation as an exterior differential 
system. The simplest example, but by no means the only example, is the wave 
equation u xy = 0. The wave equation can be encoded as adifferential system on 
manifold of dimensions 5, 6 and 7 by 

A v 

X\ = { 9, dp A dx, dq A dy }, V\ — { 9, dp, dx}, V = { 9, dq, dy }, 

! 2 ={0,9 x ,dqAdy}, V 2 = {9,9 x ,dp,dx}, V 2 = {9,9 x ,dq,dy}, 
1 3 = {9, 9 x ,9 y }, V a = { 9, 9 X , 6 y , dp, dx },V 3 ={ 9, 9 x ,9 y , dq, dy }, 

where 9 = du — pdx — q dy, 9 X = dp — r dx and 9 y — dq — t dy. Each of these 

satisfy the definition of a decomposable, Darboux integrable differential system. 

. . a ^ 

Consistent with Corollarv l2.6i one easily checks that li — Vi FR W 

We remark that for the standard contact system on J 1 (R, R), that is, 

Xt = {9, dp A dx + dq A dy}, the Pfaffian systems V4 — { 9, dp, dx } and V = 

{ 9, dq, dy }, are singular systems which form a Darboux pair but V ^ FR V4 7^ I4. 

Equality fails because I4 is not decomposable. | 

Remark 2.12. In very special cases, such as Liouville's equation u xy = e u , 
one finds that V C\ V = V nV^ = {0}so that the vector space sums 
appearing in (II. 8|) are direct sum decompositions. Under these circumstances, 
our 0-adapted coframe is simply {6,&, &}. This initial coframe is then auto- 
matically 2-adapted (see Theorem l4.3[) and the sequence of coframe adaptations 
needed to prove Theorem 14 . 1 1 can begin with the third coframe adaptation (Sec- 
tion 4.2). 

Remark 2.13. We define an involution of a Darboux pair { V,V} on M 

O A. V 

to be a diffeomorphism $: M — > M such that $ 2 = idjvf and <f>*(V r ) = V. For 
such maps $*(V K ') = V and therefore, given 1 adapted co-frame elements 
fj and <x, one may take 

17 = $*{fj) and & = $*(<t). 

Such involutions are quite common, especially for differential equations arising 
in geometric and physical applications, where the diffeomorphism $ arises from 
a symmetry of the differential equations in the independent variables. 
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3 The group theoretic construction of Darboux 
pairs 

3.1 Direct Sums of Pfaffian Systems and Trivial Darboux 
Pairs 

Let Wi and W2 be differential systems on manifolds Mi and M 2 . Then the 
direct sum Wi + W2 is the EDS on Mi x M2 generated by the pullbacks of Wi 
and W2 by the canonical projection maps m : M\ x M2 — > Mi. The following 
theorem shows that if Wi and W2 are Pfaffian, then Wi + W2 is Darboux 
integrable in the sense of Definition 11.31 

Theorem 3.1. Let W\ and W2 be Pfaffian systems on manifolds Mi and M2 
such that = W 2 (oo) = 0. Then the direct sum J = Wi +W 2 on M x x M 2 

is a decomposable Pfaffian system. The associated singular Pfaffian systems 

W = Wi® A 1 (M 2 ) and W = A 1 (Mi) © W 2 (3.1) 

form a Darboux pair and 

Wi + W 2 = W ffl VV, (3.2) 
where Wi, W2 5 VV, W are i/ie differential systems generated by W\, W2, W, W . 

Proof. It is easy to check that J is decomposable with singular systems (|3.ip . 
Properties [i] and [ii] of Definition 1 1.21 follow directly from l|3.ip and the simple 
observations that 

W (ca) = A 1 (M 2 ) and W (oo) = A 1 (M 1 ). (3.3) 

Let uj — Yliifi ®\ + 9i ^2); where the coefficients fi, gi G C°°(Mi x M 2 ) and 
the 9\ € are 1-forms in J. Then to check property [iii], it suffices to note 
that every summand in the exterior derivative 

doj = ^d 1 f l A e\ + d 2 f l a e\ + ^ d x e\ + d l9j a e J 2 + d 29] a e{ + 9l d 2 e l 2 , 

i 

where d,; = dM t , lies either in Cl 2 (W) or Q 2 (W). Equation (|3.2p is a direct 
consequence of Theorem 12.61 | 

Darboux pairs of the form (|3.1|) are said to be trivial. 
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3.2 Symmetry Reduction of Darboux Pairs 

A Lie group G is a regular symmetry group of a differential system 

X on a manifold M if there is a regular action of G on M which preserves X. By 
the definition of regularity, the quotient space M /G of M by the orbits of G has 
a smooth manifold structure such that [i] the projection map q: M — * M/G is 
a smooth submersion, and [ii] local, smooth cross-sections to q exist on some 
open neighborhood of each point of M/G. The G reduction of X or the 
quotient of X by G is the differential system on M = M/G given by 

I/G-{weO*(M)|q*(w)eI}. (3.4) 

Details of this construction and basis facts regarding the G reduction of differ- 
ential systems are given in [Tj. 

To obtain the main results of this section, we shall require that G be trans- 
verse to X in the sense which we now define. Let I be any G invariant sub- 
bundle of T*M and let V be the infinitesimal generators for G. We say that G 
is transverse to I if for each x £ M 

r.nannJ, = {0}. (3.5) 

We say that G is transverse to the differential system X if G is transverse to the 
sub-bundle of T*M whose sections are 1-forms in X. With this transversality 
condition, X/G is assured to be a constant rank differential system whose local 
sections can be determined as follows. If A is any sub-bundle of A k (M) or sub- 
algebra of A*(M), then the sub-bundle or sub-algebra of semi-basic forms in 
A is defined as 

A sh = {uo e A\ X-' uj = for all X e T}. (3.6) 

Denote by 2" s b the sections of X which are semi-basic, that is, which take values 
in A* b (M). Then one can show (see Lemma 4.1 of [T]) that 

(I/G)p=t*(I ub \ u ), (3.7) 

where U is any (sufficiently small) open set of M , £ : U — > M is a smooth 
cross-section of q, and U = q _1 (?7). We remark that the G reduction X/G of a 
Pfaffian system X need not be Pfaffian (but see [c] below) . 

For the proof of the next theorem on the G reduction of Darboux pairs, we 
shall require the following elementary facts regarding the reduction of Pfaffian 
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systems. Here / and J are Pfaffian systems with regular symmetry group G, we 
assume that G is transverse to /, and the G bundle reduction of / is defined as 

I /G = {uj e A 1 (M) | q*(uj) G /}• (3.8) 

By transversality, this is a sub-bundle of T* M of dimension dim / — q, where q 
is the dimension of the orbits of G. 

[a] If / C J, then G is transverse to J. 

[b] If G is transverse to 1 J, then (I/G) (J/G) = (I n J)/G. 

[c] If G is transverse to the derived system then the quotient differen- 
tial system I/G is the constant rank Pfaffian system determined by I/G and 
(I/G)' = I'/G . 

[d] If G is transverse to jM, then (I/G)^ = I^/G. 

Facts [a] and [b] are trivial, while [c] is Theorem 5.1 in pQ. Fact [d] follows 
from [a] and [c] by induction. 

Finally, we shall also use the following technical observation. If Z is a sub- 
bundle of/, then the canonical pairing (X, lu) = uj(X) on T x x Z x , where x € M, 
is non-degenerate if and only if 

I = cf{I/G)®Z. (3.9) 

Such bundles Z can always be constructed locally on G invariant open subsets 
of M. 

A v 

Theorem 3.2. Let { VF, W} be a Darboux pair on M. Suppose that G is a 

A v 

regular symmetry group of W and W and that, in addition, G is transverse to 
WnW (oo) and W (oo) nW. Then the quotient differential systems { W/G, W/G } 
on M are (constant rank) Pfaffian systems and form a Darboux pair. 

Proof of Theorem^ Since W {Qo) nW C W {oo) C W' C W and G is transverse 
to W {co) n W, it follows that G is transverse to W {oo) , W' and W. Thus, by 
[c], V = W/G and, similarly, V = W/G arc the Pfaffian systems for V = W/G 
and V = W/G. By [d] 

(W/G) (oo) = W {oo) /G and [W/G) {oo) = W (oo) /G. (3.10) 
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We now verify the three conditions in Definition 11.21 for { V, V } to be a 
Darboux pair. From (|3 . 8[) . (|3 . 10[) and fact [b] we deduce that 

A A 

dim(F) = dim W — q : 
dimO> (oo) ) = dim( W {oo) /G) = dim W (oc) - q, and 
dim(y n V (oo) ) = dim ((W/G) n {W (oo) /G)) = dim ((W fl VF (oo) )/ G ) 
= dim(VFrW (oo) )-g. 

From these equations and property [i] in the definition of Darboux pair (applied 

A 

to { W, W }) we calculate 

dim (V + V (oo) ) = dim(F) + dim(F (oo) ) - dim(F n F (oo) ) 

= dim W + dim VF (oo) + dim(VFnlF (oo) ) -9 
= dim M — q = dim A/. 

This proves that V + V {oo) = T*M. Similar arguments yield V {oo) + V = T*M 
and so {V V} satisfy property [i] of a Darboux pair. 

A v 

To check that {V, V} satisfy property [ii] of a Darboux pair, we use [b] and 
(|3.10p to calculate 

n ^(00) = ( ^(-) /G) n iw [co) /G) = (w [co) n w [co) )/G = 0. (3.11) 

To prove property [iii] of a Darboux pair, we first prove that 

[A 2 (VF) + A 2 (W)] sh = A 2 (W sh ) + A 2 (W sb ). (3.12) 

A (qqS V 

It suffices to check this locally. Since G acts transversally to W fl W, there 
is a locally defined sub-bundle Z of W n VF such that the canonical pairing 
on r x Z is pointwise non-degenerate. For any such sub-bundle Z we have that 
W = W sb © Z and W = W sb Z and therefore 

A 2 (VF) + A 2 (FT) = A 2 (FF sb ) + A 2 (FF sb ) + T*M A Z. 

This (local) decomposition immediately leads to (|3.12|) . To complete the proof of 
property [iii], let ui G FnF. Then q*(w) G VFlH VF and hence, by the properties 
of the Darboux pair { W, W} and (IXT21) . q*(dw) G ^ 2 (TT sb ) + ^ 2 (VF sb ). We 
pullback this last equation by a local cross-section of q and invoke (|3.7[) to 
conclude that duj G fi 2 (F) + VL 2 (V). | 
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Theorem 3.3. Let { W, W} be a Darboux pair on M. Suppose that G is a 
regular symmetry group of W and W and that, in addition, G is transverse to 
W ( °° ] nW andWn W (oo) . Then 

(VV ffl VV)/G = (VV/G) ffl (VV/G). (3.13) 

Proof. We first remark that Theorem E2] implies that V = W/G and V = W/G 
define a Darboux pair and, second, that it suffices to check (]3 . 13[) locally. 

Transversality implies that there are locally defined sub- bundles Z C W n 
W and Z cWn W (oo) on which the natural pairing with T (the infinitesimal 
generators for G) is non-degenerate at each point x £ M . It follows from (|3.9p 
that 

W ioo) nW = q*(V (cc) nV)®Z, # (oo) = q*(V {oo) ) © Z, (3.14a) 

#rw (oo) = q *(yn^ (oo) )©z, w> (oo) = q *(^ (oo) )ez, (3.14b) 

# = q*(y)©Z, VV = q*(1>)eZ, and (3.14c) 

VFnw = q*(7nf) + z = q*(yn£) + z. (3.i4d) 

If £ and £ are local basis of sections for Z and Z then, because £ e It and 
C S W { °°\ dC = mod VV (oo) and d£ = mod VF (oo) , or 

<iC = mod {q*(f> (oo) ), C} and dC, = mod { q*(V (oo) ), C }. 

These equations, together with (|3.14c|) , show that the Pfaffian systems W and 
VV are algebraically generated as 

VV = {q*(V), <} and VV = { q*(V), C}. (3.15) 

To complete the proof of the theorem, we now make the key observation that 
since £ g W n W , (|3. 14c|) implies that we may express these forms as linear 
combinations of the £ and the forms in q'(^ny). Hence VV is also algebraically 
generated as VV = {q*(V),£}. We may list an explicit set of local generators 

A v 

for q*(V) and q*(V) using the 1-adapted coframe constructed for the Darboux 
pair { V, V } in Theorem 12.91 The arguments given in the proof of Theorem 
12.101 can then be repeated to show that VV ffl VV = { q*(V ffl V), C } and hence 
(VV FR VV) s b = q*(V ffl V), as required. | 
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Corollary 3.4. Let W\ and W% be Pfaffian systems on manifolds M% and Mi 
with w[°° ] = W 2 (oo] = 0. Suppose that 

[i] a Lie group G acts freely on M\ and Mi and as symmetry groups of both W\ 
and Wi; 

[ii] G is transverse to both W\ and W%; and 

[iii] the diagonal action of G on M\ x Mi is regular. 
Then the quotient differential system 



is a decomposable, Darboux integrable, differential system with respect to the 
Darboux pair 



Proof. As in Theorem 13.11 define the Darboux pair W — W\ © A 1 (Mi) and 
W = A 1 (M 1 ) W 2 . To conclude from Theorem [O that U = W/G, U = W/G 
are also a Darboux pair, we observe that [i] and [ii] imply that the diagonal 
action of G on Mi x M 2 is transverse to the Pfaffian systems W\ and Wi pulled 
back to Mi x M 2 (a fact which is not true without [i]. We also note that, by 



J = (Wi+W 2 )/G 



(3.16) 



U = (Wi © A 1 (M 2 ))/G and U = (A 1 {M 1 ) © W 2 )/G. 



(3.17) 



(Els , 



WnW 



Wi and W' 



(oo) 



n w = w 2 . 



(3.18) 



We now use Q3.2[) and Theorem 13.31 to deduce that 



J = (Wi + Wi)/G = (VV ffl W)/G = (W/G) ffl (W/G) =UmU. 



An application of Corollary 12 . 101 completes the proof. 
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4 The coframe adaptations for a Darboux pair 

In this section we present a series of coframe adaptations for any Darboux 

A v 

pair { V, V}. These coframe adaptations lead to the proof of the following 
theorem. 

Theorem 4.1. Let {V, V} be a Darboux pair of Pfaffian differential systems 
on a manifold M. Then around each point of M there are 2 coframes 

{6\n a ,fr a } and {8 l ,n a ,Tc a } (4.1) 

such that 

V = {8\n a }, V°° = {n a }, V = {9\n a }, V°° = { n a } (4.2) 

and with structure equations 

dff 1 = -G^gX" A n 13 + ^C' ,0 J A § k and 

? i (4-3) 

d9 i = -H i ab Ti«At b --C i jk VA9 k . 

The coefficients C* k are the structure constants of a Lie algebra whose isomor- 

A V 

phism class is an invariant of the Darboux pair { V, V }. 

We have already constructed the first coframe adaptation in Section 2. The 
second coframe (Section 4.1) adaptation provides us with a stronger form for 
the decomposition of the structure equations than that initially given by (|2.5p . 
The third and fourth coframe adaptations (Section 4.2 and 4.3) lead to the 

A *v- 

definition of the Vessiot Lie algebra Qess(V, V), with structure constants Cj fe , 
as a fundamental invariant for any Darboux integrable differential system. The 
importance of this algebra was first observed by Vassiliou who introduced 
the notion of the tangential symmetry algebra for certain special classes of 
Darboux integrable systems. The tangential symmetry algebra is a Lie algebra 
of vector fields T on M which is isomorphic (as an abstract Lie algebra) to the 
Vessiot Lie algebra. The tangential symmetry algebra also plays a significant 
role in Eendebak's projection method [TD] for Darboux integrable equations. 

Of course, the tangential symmetry algebra T determines a local transforma- 
tion group on M but this transformation group will generally not be the correct 
one required for constructing the superposition formula (See Example 6.1). One 
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final coframe adaptation (Section 4.4) is needed to arrive at the proper trans- 
formation groups (See Definition 15. 7|) for this construction - this is the coframe 
that satisfies the structure equations given in Theorem 14. 11 

The coframe adaptations 1-4 involve only differentiations and elementary 
linear algebra operations. For the fifth coframe adaptation, in the case where 
the Vessiot algebra is not semi-simple, one must integrate a linear system of 
total differential equations (see (|4.61[) ) and also solve the exterior derivative 
equation d(a) — (3, where (3 is a closed 1 or 2-form (see (|4.64|) and (|4.73|) ). 

4.1 The second adapted coframe for a Darboux pair. 

A v 

Let {V , V} be a Darboux pair and let {9, <r, fj, <r, 17} be a 1-adapted 
coframe with dual frame {d s , d$., d^, d„. d^}- From this point forward, we 
are solely interested in adjustments to the forms 6 that will simplify the 1- 
adapted structure equations 

dO = Afr Act + B a Act mod { 0, fj, fj }. (4.4) 

Of the "mixed" wedge products 

ctAct, <t A 17, r)A(T, and fj A 77 (4.5) 

the products & A <x are the only ones definitely not present in (|4.4|) . We shall 
now show that it is possible to make an adjustment to the 1-forms 8 of the type 

9' = 8 + Pfj + Qf) (4.6) 

so that the structure equations for the modified forms 0' are free of all the wedge 
products (|4.5[) . We begin with the following simple observation. 

Lemma 4.2. //{ d g , d~, d~, d„. } is the dual frame to a 1-adapted coframe 
for the Darboux pair { V, V }, then 

[d~,d^}=0. (4.7) 

Proof. It suffices to note that none of the structure equations for a 1-adapted 
coframe contain any of the wedge products <x A <r. | 

The construction of the second adapted coframe is completely algebraic in 
that only differentiations and linear algebraic operations are involved. Let S a = 
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<9-o and S a = da a and define two sequences of vector fields inductively by 

A A A ■■■ f 

S aici2---ai — [ ^aiaz-'-ag— 1 i ^a?: ] and S aia2 ... a ^ — [ S aia2 .. .q^ ^ , S a ^ ] . (4.8) 

Because the vector fields S ae G , the vector fields Sa^. .. af belong to H and 
are therefore linear combinations of the vector fields {9 e , 9^., 9~}. However, 
because da = 0, one can deduce by a straightforward induction that in fact 

S ai a 2 ---a e € span{9 e , d~ } and, likewise, S^aa-a* S span{9 e , d^}. (4.9) 

The Jacobi identity also shows that 

ff (oo) = span{ S aia2 ... ae } e >! and if (oo) = span{ S ai0t2 ... ae }e>i. (4.10) 

From (|4.7p it also follows that 

A ^ 

[ S aia2 ... ak , *S , Q 1 Q 2 ---a^ ] = 0* (4-H) 

for k,l > 1. By virtue of (|4.10p we can choose a specific collection of iter- 
ated Lie brackets S ai a 2 ---a e and S , ctl a 2 ... Qf , t > 2, which complete {9 e , 9^. } 
and {9 e , c?^. } to local bases for i/*- 00 -* and j/^ 00 -* respectively. Denote the 
iterated brackets so chosen by {§'} and {§'}. Since span{9 e , d~, 9~} = 
span{9 e , d^, §'} and span{9 e , d^, 9^} = span{9 e , 9^, s'}, this gives us a 
preferred 0-adapted frame { d e , 9- , s , , s' } where, on account of (|4.11[) . 

[9-, 5 ] = O, [O fr ,« / ] = [s',dz} = 0, [s',s']=0. (4.12) 

It follows that the coframe { 9', <r', f\ , <x', 17' } dual to { 9 e , 9^., s', 9,j, s' } is 
a 0-adapted coframe for the Darboux pair {V, V}. Equations (|4.9|) imply (14. 6|) . 

On account of (|4.12[) . the structure equations for the forms { 0' } are free of 
all the wedge products a A &' , &' Af]', f) A a', and f)' hi)' . We can express 
this result by writing 

de'e n 2 (v {QC) ) + n 2 (v {QC) ) mod{0'}. 

The coframe { 9', <x, fj, <r, f\ } therefore satisfies the criteria of the following 
theorem. 

A *v* 

Theorem 4.3. Let { V, V} be a Darboux pair on a manifold M. Then about 
each point of M there exists a 1-adapted coframe { 6, <r, r), <r, 17 } with structure 
equations (I2.13P and 

d9 = Atv Att + JStt Att mod{0}, (4.13) 

where tt and -ft denote the tuples of forms (o - , 17) and (er, fy). 
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Coframes which satisfy the conditions of Theorem 12.81 are call 2-adapted 
coframes. 
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4.2 The third adapted coframes for a Darboux pair 

Written out in full, the structure equations (|4. 13|) are 

d9 = A7rA7r + BTTATT + COAO + M7i-A0 + NTTAO. (4. 14) 

We obtain the third coframe reduction for a Darboux pair by showing that a 
change of coframe 9' = P 6 can be made so as to eliminate either all the wedge 
products 7T A 6 or all the wedge products tt A in (|4.14p . The construction of 
this coframe uses another set of iterated Lie brackets. 

We start with a 2-adapted coframe {8, <x, fj, er, fj,}, where 7r = ( &, fj ) and 
7r = (<t, 77), and introduce the provisional frame {6, 1,1}, where L = dl and 
1= dl. This is not a 0-adapted coframe although it is the case that 

span{ t } = span{ tv } and span{ I } = span{ it }. (4-15) 

From (14. 14[) and (|4.15|) we find that the structure equations for this coframe are 

dt = 0, dl = 0, and 

(4.16) 

dG = AiAl + BLAL + C6A0 + MtA6 + NiA6. 
Denote the dual to this provisional coframe by { d e . U, U}. Since 

span{ d~ } = axmV — ann{0, t, <x} C ann{#, 1} = span{ 17}, 
we have that 

dg. C span{ U } and C span{ U }. (4-17) 

From (|4.16p it follows that the structure equations for the vectors fields U 
and U are 

[U, U]=-Adg, [U,d g ]=-Md , 

[TJ,tJ}=-Bd g , [U,d e ]=-Nd e , (4.18) 
and [U,U]=0. 
As in Section 4.2, define two sequences of vector fields inductively by 

A A A v A A 

U aid2 •••0,1 — [U aia.2 - ■ -a,£_i : U a ^ ] and U ai a 2 — [ U ai0t2 ■■■a ( ,_ l U Q£ ]. (4.19) 

A simple induction argument, based upon the last of (|4.18[) and the Jacobi 
identity, shows that 

[U ai a 2 ... ak i t^Q 1 Q 2 '"CK£ ] — ^ - (4.20) 
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for k, I > 1. On account of (f2~T|) , |2.9p and (|4.17|) the iterated brackets U aia2 ... at 
will span all of H and therefore, by the first two equations in (|2.9p , we can choose 
a basis X for X = H^ 00 ^ n iy' 00 -* (see (I2.7p ) consisting of a linear independent 
set of the vector fields U ai a 2 ---a e - Alternatively, we can choose a basis Y for 
K consisting of a linear independent set of the vector fields U ai a 2 ---a t - Due to 
(|4.20p these two bases for K satisfy 

[X,U}=0, [X,Y]=0, [Y,U] = 0. (4.21) 

Denote the dual of the coframe {X, U, U } by { 0x, t , I } and the dual of the 
coframe {Y, U, U} by { 0y, £ , £ }. 

/\ ^ 

Theorem 4.4. Let { V, V} be a Darboux pair on a manifold M. Then about 
each point of M there are two 2- adapted coframes 

[6 X , <r, fj, ct, fj} and {0 Y , <r, fj, a, fj} 

with structure equations (|2.13p . 

d6 x = AtvAtv + BtcAtc + C0xA0x + MtvA0x, (4.22) 

and 

d0 Y = ETvATv + FnAn + K0 Y A0 Y + N^A0 Y . (4.23) 

Moreover, span{$x} = span{0y} and the vector fields X, Y belonging to the 
dual frames { X, d-, 9-, } and { Y, d- , d~, d% , 5^ } safe/?/ 

[X,Y}=0. (4.24) 

Coframes which satisfy the conditions of Theorem 14.41 are called 3-adapted 
coframes. 

Remark 4.5. In the special case where the pair { V, V} admits an involution 
(see Remark 12. 13[) the forms Y may be defined by Y — $*(0x)- Then the 
structure equations (|4.23p can be immediately inferred from (|4.22p . 
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4.3 The fourth adapted coframe and the Vessiot algebra 
for a Darboux pair. 

We now show that any 3-adapted coframe may be adjusted so that the 
structure functions Cj fc and Kj k (see (|4.22[) - (|4.24p ) are constants and satisfy 

r ii — — jc i 
u jk - » ■ 

A V 

Theorem 4.6. Let {V, V} be a Darboux pair on M . 

[i] Then, about each point of M , there exist local coframes { 6 X , TT a , rr a } and 
{9\r, : K a ,TT a } which are 3-adapted and with structure equations 

dO l x = \K b tt- A n b + 7r Q A 7r^ + ic7j,4 A 0| + M^ a A 4 (4.25) 

and 

dOV = \K b 7r° A 7r fc + l -Fl n<* A ^ - A # + A^tP A 4- (4-26) 

T/ie structure functions Cj fe = — • are constant on M and are the structure 
constants of a real Lie algebra. 

[ii] The isomorphism class of the Lie algebra defined by the structure constants 
C* fc is an invariant of the Darboux pair { V, V}. 

Coframes which satisfy the conditions of Theorem l4.6l are called ^-adapted 
coframes. As with the 2 and 3-adapted coframes defined previously, 7r = 
((T,fj) and 7r = (<r,i7), with structure equations (|2.13[) . In passing from the 
3-adapted coframes to the 4-adapted coframes only the forms 9 X and B % Y are 
modified. The (abstract) Lie algebra defined by the structure constants Cj fe 
is called the Vessiot algebra for the Darboux pair { V, V } and is denoted 
by Oess(V, V). The vectors fields Xi in the dual basis { Xj, U a , U a } to the 4- 
adapted coframe { , 7r a , n a } define a realization for X>css(V, V), with structure 
equations [Xj,X k ] = -C] fe X 4 . 

Proof of Theorem \4-6\ . Theorem 14.41 provides us with two locally defined 3- 
adapted frames { Xj, C/ a , [/ a } and { Yi, U a , U a } and dual coframes { 6 X , n a , 7r Q } 
and {6>^,7T a ,7f Q }. The structure equations are (j4~22]) — (|4~24]) . We begin the 
proof of part [i] by showing that the Cj fe are functions only of the first integrals 
I a while the K l - k are functions only of the first integrals I a . 

Since span{ 6x} = span{ 6y}, there is an invertible matrix Q such that 

6 X = Q[>){ . (4.27) 
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On account of the structure equations (|4.22p and the fact that the vector fields 
Xi and Yj commute, the identity 

d# x (X jt Y k ) = - - O x ([X j; Y k }) 

leads to 

XiiQi) (4.28) 

We shall use this result repeatedly in what follows - it is equivalent to the fact 
that the vector fields X{ and Yj commute, a property of the two 3-adapted 
coframes which is not encoded in the structure equations (|4.22[) and (14.23(1 . 

We next substitute (|4.27[) into (|4.22j) and equate the coefficients of 9 3 Y A fly- 
to deduce that 

XtiQDQt - XdQ^Qi + Q\K\ k = C\ m Q\Q^. 

By virtue of (|4.28[) . this equation simplifies to 

Q\K) k = -C\ m Q\Q%. (4.29) 

Also, by equating to zero the coefficients of tt " A 9 X A 9 X and 7r° A fly A fly in 
the expansions of the identities d 2 x = and d 2 l Y = 0, wc find that 

U a {C) k ) = Q and U a (K} k ) = 0. 

By Lemma 12.81 these equations imply that 

C) k e Int(V) and K) k € Int(F). (4.30) 

Our next goal is to show that the coframes Q l x and fly may be adjusted so that 
Kj k = —Cj k while still preserving (|4.22[) — (I4.24[) . To this end we introduce local 
coordinates (I ,I a , z rn ) on M satisfying (I2.10p and (12. 1 1[) . Then, on account of 
(|4.30[> . equation (|4.29|) becomes 

Qi(I a J a ,z m )K*(I a ) = -^(/ a )Q^(/ Q ,/ a ,z m )Q J fc (/ a ,/ Q ,z m ). (4.31) 

Evaluate this equation, first at a fixed point (/q,/q, z™) and then at the point 
(7 a ,^,z m ). With 

K%=K*{TZ), C% = C^{ll), Ql=Ql(I a ,I a ,z™), and 

Ql(I a ) = Ql(I a J« lZ ™) 
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the results are 

QiK% = -CL QiQj and Qi(I a ) k% = -Ci k {I a ) Q*{I a ) Q]{I a ). (4.32) 
It then readily follows that the matrix 

Pi(I a ) = Q}(I a )(Q-% (4.33) 

satisfies 

HP)C% = Ci k {I a )Phl a )Pf{I a )- (4.34) 
The 1-forms Q x defined by 

o x p; = e x (4.35) 

then satisfy structure equations of the required form (|4.25[) , where the structure 
functions Cj fc coincide with the constants C) k - Finally, by evaluating (|4.31[) at 
(Jq, I a , z™), it follows that 

R){I a ) = (Q-\llT,zZ))) (4.36) 

satisfies 

-Rl{I a )C% = Kf lk (I a ) R'l(J a ) R*(I a ) (4.37) 
and the 1-forms Q % Y defined by 

y> Y R\ = B Y (4.38) 

satisfy the required structure equations (|4.26p . again with C l - k = C) k - 

The next step in the proof of [i] requires us to check that the coframes we 
have just constructed are still 3-adapted. We must therefore show that the dual 
vector fields Xi and Yj commute and for this it suffices to show, because of our 
remarks at the beginning of this proof, that 

X^Ql) = C l 3l Q{ where X = Q)6 Y , (4.39) 

The proof of this formula follows from three simple observations. We first 
note that (|4351) and imply Q\ = (P-^Q^R^. Secondly, because 

{ X 1 , % a , n a } is the dual frame to { 6 X , % a , n a }, we have Xi = P^Xj. Finally, 
because P] = P*(I a ) and R) = R){I a ), it follows that Xj{P l k ) = Xj{R k ) = 0. A 
straightforward calculation based upon these three observations and equations 
P~2"5|) and P~M]) leads to (HTM)) , as required. 
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Finally, by equating to zero the coefficients of 6 X A (P x A 9 X in the expansion 
of the equations d 2 9 x — one finds that the Cj k satisfy the Jacobi identity and 
are therefore the structure constants of a real r-dimensional Lie algebra. 

To prove [ii], let { 9 X , Tt' a , 7r' Q } be another coframe adapted to the Darboux 
pair { V, V } with structure equations 

d9 x = V a A n' b + ifl*, Tr' a A n" 3 + -C^ A 9 ,k + M^)n' a A Q' j '. (4.40) 

From the definition of a 0- adapted coframe (see (|2.6bj) ) we have 

e x = r)& + s l p ?) p + s l q 7f, 

where the matrix Rj is invertible. Substitute this equation into (|4.40p and then 
substitute from (|2.13| and (|4.25| . From the coefficients of & A 6 we deduce that 
d~Rj = 0. By Lemma I2T51 this implies that Rj S Inv(V') in which case one finds 
from the coefficients of 9^ A 9 k that 

Cli r>l nm r>if~tl 
lm n j n k — n l^jk- 

Hence the structure constants Cj fc and C'j k define the same abstract Lie algebra 
and the proof of part [ii] is complete. | 

Corollary 4.7. Let {V , V} and {W,W} be Darboux pairs on manifolds M 
and N respectively and suppose that <f>: M — > N is smooth map satisfying 

<j)*{W)<zV and <j)*{W)dV. (4.41) 

Then <fi induces a Lie algebra homomorphism 

4>: otSs(V,V) -> vess{W,W). (4.42) 

Proof. The proof of this corollary is similar to that of part [ii] of Theorem 14.61 

A v 

Let { 9 x ,ir' c , tt' 1 } be a 4-adapted coframe on N for the Darboux pair { W, W }, 
with structure equations 

d9 x = -A% n' c A n' d + ^Bfi A n' 5 + ^K s rt 9' r A 9 n + M£n' a A 9 n . (4.43) 

The constants K' t are the structure constants for the Lie algebra oess(W, W). 
Denote the dual frame by { X' r , U' c , U'^ }. 

Let { 9 X , 7r a , 7r Q } be a 4-adapted coframe on M for the Darboux pair{ V, V } 
with the structure equations (|4.25[) . The inclusions (|4.41[) imply that 

cf>*(Wn W) C V(1V, (j)*(W ico) ) dV {oo) and (j)*(W {oo) ) C v {oo) 
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and therefore there are functions Kf, S p , S q T a and T a on M such that 

r(0'i) = Rtei c +S s p fr + S s q fj«, ^(^' c )=f c a n a and cj>*{^) = f 

We pullback ()4.43p using these equations and substitute from (|4.25|) . The same 
arguments as in the proof of [ii] of Theorem 14.61 now show that Rf 6 Inv(V) 
and RiCj k — Ktf.RjR\.. This proves that the Jacobian mapping (f>*(Xi) = R\X' t 
induces a Lie algebra homomorphism of Vessiot algebras. | 

Remark 4.8. Let {0 X , ^ a } and { 8y , 7r a , tt a } be two coframes which are 
4-adapted to the Darboux pair { V, V } and which satisfy the structure equations 
(|4.25[) and (|4.26[) . Then it is not difficult to check that the commutativity of 
the dual vector fields Xi and Yj is equivalent to the supposition that the change 
of frame (|4.27[) defines an automorphism of the Vessiot algebra, that is, 

Q'i^ k = C\ m Q\Q^. (4.44) 

In the next section we shall need all the derivatives of the matrix Qj . The w a 
and 7r a components of dQ^ are easily determined by substituting (|4.27p into 
(|4.25| and comparing the result with (|4.26|) . If we then take into account (|4.28|) 
and (|4~44]) we find that 

dQ) = Q\M\ - Q\Nj + q,Qi4, (4.45) 

where Mj = M % a ^k a and iVj = N^jfr a . We note, also for future use, that 

Ai b = Q)Ei b and B^ = Q)F^. (4.46) 

where E 3 ab and are defined by (12. 13|) . 
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4.4 The fifth adapted coframe for a Darboux pair and the 
proof of Theorem 14.11 

Let { V, V } be a Darboux pair on M and let { 6 X , n a , n a } and { 9 l Y , n a , n a } 
be local coframes on M which are 4-adapted to the Darboux pair { V, V } and 
which therefore satisfy the structure equations (|4.25[) and (|4.26[) . In this section 
we shall prove that it is possible to define forms 

9 l = k)6 x + fa n a and ¥ = k)e Y + i>l n a , (4.47) 

where R), fa <E Int(F) and R], ^ € Int(V), 
which satisfy structure equations 

d6 l = -G^ w a A ttP + -C}J j A 9 k (4.48) 

and 

dO l = \K b 7T a A l -C) k fr A 9 k . (4.49) 

These are the structure equations for the Darboux pair {V, V} announced in 
Theorem 14.11 

We shall focus on (|4.48p and simply note that the proof of (14.49)) is similar. 
Our starting point are the equations (|4.25|) and (|4. 26[) . in particular, 

dO x = \K b * a A * b + \K& tt q A ^ + \c) k Q x A 6 X + M^ a A 9 X . (4.50) 

In what follows it will be useful to introduce the 1-forms and 2-forms 

M] = M e a] n a , N? = A^-tP, A 1 = A l ab % a A £ b F l = F^rc a A % . (4.51) 

By setting to zero the coefficients of 7r a A 7r Q A 9 X and 7r Q A 7r a A 7r b in 
the equations d 2 8 x = 0, we find that U a (AP aj ) = and U a {A l ah ) = 0. These 
equations imply that M^j € Int(V) and A l ab 6 Int(F) and therefore (see Lemma 
I2~5D 

dM* aj =U b (AP aj )n b and dA\ h = C/ c (4 6 ) tt c . (4.52) 

Bearing these two results in mind, we then respectively deduce from the coeffi- 
cients of n a A 6 X A 6 X , n a A n b A J X and n c A n a A 7T 6 in d 2 9 x = that 

M^Cj fc + M f ak q e - = 0, (4.53) 

dM) - M\ A Mj + C^A* = 0, and (4.54) 
dA l - A l M\ = 0. (4.55) 
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The proof of Theorem 14.11 hinges upon a detailed analysis of equations (|4. 53[) 
(|4.55|) . We deal first with (|4.53[) since this is a purely algebraic constraint. It 
states that for each fixed value of a, the linear transformation M a : Oess — > oess 
defined by 

M a (X.j) = M^Xi 

is a derivation or infinitesimal automorphism of the Vessiot Lie algebra ttess = 
Oess(y, V). Consequently, to analyze this equation we shall invoke some basic 
structure theory for Lie algebras. Specifically, we shall consider separately the 
cases where Oess is semi-simple, where Oess is abelian, and where Oess is solvable. 
Then we shall make use of the Levi decomposition of Oess to solve the general 
case. 

Case I. We first consider the case where the Lie algebra Oess is semi-simple. 
Here the proof of (|4.48j) is relatively straightforward and is based upon the fact 
that every derivation of a semi-simple Lie algebra is an inner derivation (see, for 
example, Varadaraian[2"5]. page 215). In fact, because the proof of this result is 
constructive, we can deduce from (|4.53[) that there are uniquely defined smooth 
functions 4> a € Int(V) such that 

M^=tf a C\ y (4.56) 

The forms Q % defined by 

& = #x + ^ a (4.57) 
then satisfy structure equations of the form 

dP = \jV ab n a A n b + \w a(i n a A ^ + \c) k fr A 9 k . (4.58) 

For these structure equations, equations (|4.54|) reduce to C\^A l ah = 0. Since we 
are assuming that oess is semi-simple, the center of Oess is trivial and therefore 
this equation implies that A l ab — 0. The structure equations (|4.58j) then reduce 
to the form (|4. 50[) . as desired. 

Case II. Now we consider the other extreme case, namely, the case where 
Oess is abelian. Strictly speaking, we need not treat this as a separate case but 
the analysis here will simplify our subsequent discussion of the case where Oess 
is solvable. When Oess is abelian the structure equations (|4.50[) are 

d0\ = \A l ah n a A n b + ^ A ^ + M l aj n a A ^, (4.59) 
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equation (|4.53[) is vacuous and, by virtue of (|4.52p , (14.54[) simplifies to 

dM) - M\ A Mj = 0, where M) = AP aj n a . (4.60) 

By the Frobenius theorem we conclude that there are locally defined functions 
R) G Int(V') such that 

d{R)) + R\ M] = and det(R)) ^ 0. (4.61) 

(For this precise application of the Frobenius theorem see, for example, Flanders 
[TTj . page 102.) We remark that this is the first (and only) step in all our coframe 
adaptations that require the solution to systems of linear ordinary differential 
equations 

The forms 9 l = R l j0^ x are then easily seen to satisfy 

dOh = \k ab n a A n b + ^A^. (4.62) 

- - * 

Since the structure functions A\ h — R* A J ab € Int(V) it follows, either directly 

from (|4~62| or from (|4~55ll (with Mj £ = 0) that the 2-forms 

X l = \k b Z a ^Z b (4.63) 

are all d-closed. If we pick 1-forms (ff — (f^ir with S Int(F), such that 
dcjf = —X 1 ) then the forms 

e i = ei + ^ a = m j e{ + ^ a (4.64) 

will satisfy the required structure equations (|4.48[) . with C^ k = 0. By applying 
the usual homotopy formula for the de Rham complex, we see that this last 
coframe adaptation can be implemented by quadratures. 

Case III. Now we suppose that Vcss is solvable. Recall that a Lie algebra g 
is said to be p-step solvable if the derived algebras = [q < - 1 ~ 1 \ Q^ -1 '] satisfy 

0= D 0« D fl« ••• D Q^-V D fl W = {0}. 

The annihilating subspaces A^' = ann(g^) therefore satisfy 

{0} = A (0) c A (1) c A (2) • • • c A (p - X) c A (p) = g* (4.65) 

and, because the gW are all ideals, 

dA^cA^gA''). (4.66) 
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Let s = dim(g) and Sj = dimA^. If M: Q — ► Q is a derivation, then a simple 
induction shows that M: — > and therefore 

M*:AW->AW. (4.67) 

Thus the matrix representing M*, with respect to any basis for g* adapted to 
the flag (|4.65p . is block triangular. 

We apply these observations to the forms { 9 % x } and the structure equations 
(|4.50p . The forms {9 X }, restricted to the vectors Xi, define a basis for Vess* 
and consequently, by a constant coefficient change of basis, we can suppose 
that the basis { 9 l x } is adapted to the derived flag (|4.65|) . Since the notation 
for this adaptation becomes rather cumbersome in the general case of a p-step 
Lie algebra, we shall consider just the cases where Oess is a 2-step or a 3-step 
solvable Lie algebra. The construction of the coframe with structure equations 
(|4.48[) in these two special cases will make the nature of the general construction 
apparent. 

In the case where Cess is a 2-step solvable Lie algebra, we begin with a 
4-adapted coframe 

where 

span{ 9 X , . . . , 6 X } = span{ 0\, . . . , 6^ , 6 S 2 1+1 , . . . , 0| 2 } over R, 

and where { 8\, . . . , is a basis for AW(oess). In this basis the structure 
equations (|4.50p become 

dOl = ^A r ab n a At^ + ^B^F* An +MZ s Tt a A9{, and (4.68a) 

dO l 2 - \A\ h Z a A % b + n a A ^ + \C\ S 9\ l\9\ + & rj Q\ A &{ (4.68b) 

+ M l ar n a AOl + M^n" A ^. 

In these equations the indices r, s range from 1 to s% and the indices i,j range 
from si + 1 to S2' On account of (|4.66p . there are no quadratic 9 terms in 
(|4.68a|) because 9\ e A (1 \vtss) and there are no 9 l 2 A9 J 2 terms in (|4.68b|) because 
9\ e A( 2 )(Dess). There are no n a A 6\ terms in (|4.68al) by virtue of (|4.67j) . 

Since the structure equations ()4.68ap are identical in form to the structure 
equations (|4. 59[) for the abelian case, we can invoke the arguments there to 
define new forms 

9l = R r s 9{+^ a (4.69) 
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so that the structure equations (|4.68|) reduce to 

d6l = -B r afj 7f Q A tt 13 , and (4.70a) 

de\ = \k b * a A * b + \K P tP A tP + \ci s 9\ A 9{ + C rj 9\ A 9{ (4.70b) 

+ M i ar 7t a A0 r 1 +M i aj n a Aei 

It is important to track the coordinate dependencies of the coefficients in these 
structure equations as we perform these frame changes. Since the coefficients 
A\ h , Mi r and M tj in P^b]) and the coefficients R r s and <f> r a in (jCT]) are all V 



first integrals, it is easily checked that the same is true of the coefficients A z ab , 
C rs , C rj , M\ T and M\ 3 in (l4J0bl . 

The coefficients of 7r a A 7r & A 9\ and 7r a A w b A fr c in d 2 9\ — give the same 
equations as we had in the abelian case and consequently we can define 

6\ = R% + ^ a (4.71) 

so as to eliminate the 7r a A 7r b and 7r a A 9\ terms in (|4.70b[) . The structure 
equations are now 

d K = tt q A tt' 3 , and (4.72a) 

doi = \k p tt" a tt' 3 + l -c rs e\ a e[ + c), 3 e\ a P 2 + u\ r ^ a e\ . (4.72b) 

Finally, from the coefficients of 7r a A 7r b A 9\ in d 2 9\ — we find that 

rfMj = 0, where M l b = M l br n r , (4.73) 

and therefore, again by the de Rham homotopy formula, there are functions 
R r G Int(V') such that M l ar = U a {R r ). The change of frame 

9' 1 2 = 9 2 + R r 6l (4.74) 

then leads to the desired structure equations 

d$i = \Kp tt" A tt' 9 , and (4.75a) 
dh = \Kp * Q A + \c l J\ A 9\ + C rj 9[ A 9%. (4.75b) 

At this point it is a simple matter to check that the equations d 2 9\ =0 force the 

2.2. /\ 

coefficients C\ s and C r j (which belong to Int(T^)) to be constant. Equations 
(|4.75p establish (14.48)) for the case of a 2-step solvable Vessiot algebra. 
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When oess is a 3-step solvable Lie algebra we assume that the 4-adapted 
coframe { 9[, 6%, 9 Z 3 , Tr a , TT a } is adapted to the flag (|4.65|) in the sense that 

AW = span{ 6{ }, A< 2 > = span{ 61,8%} and A< 3 > = span{ 8{, 9%, 8\ }. 

The structure equations (|4.76bj) are then of the form 

d0[ = \K b tt- A n b + ~B r ap tP A **> + M r as n a A 8\, (4.76a) 

d&2 = \K b * a A * b + \Bl p n a A + \C- S 8\ A 8{ + C« 9\ A 6" (4.76b) 
+ Ml r n a A 8\ + AC * a A 9 V 2 , and 

&Q\ = \A\ h ^ + \& afi ^ (4.76c) 

+ \c\ s Bi A9{ + C l ru 9{ A 9% + 61 A 0§ + ±C*„ 9% A 2 u + C l uj 9% A 0> 
+ Ml Tt a A + M* u £ a A ^ + A/A £ Q A 4 

In these equations r, s range from 1 to s\, u, v from s\ + 1 to S2, and z, j from 
s 2 + 1 to S3. Note that the form of the structure equations (|4.76a|) — ()4.76cl) is 
preserved by changes of coframe of the type 

0i -> JSii^i + 0i ft, 2 -> R12O1 + R22O2 + and 

83 — > -Rl3^1 + R 2 302 + R33Q3 + 

.A 

where the coefficients i?y G Int(V). Exactly as in the previous case of a 2- 
step solvable algebra, we use such a change of coframe for 8\ , O2 to reduce the 
structure equations (|4.76|) to 

d6{ = \B r a ^ a h^, (4.77a) 

dO u 2 = \Kp * Q A ^ + 9\ A 8{ + C« 0* A 61 , and (4.77b) 

^^^^'A^ + ^^A^ (4.77c) 
+ \CU 61 A 6{ + Cl u 6{ A k + & rj 61 A 9> 3 + X -C l uv 6% A k + C l uj 8% A 8^ 

+ M i ar r a 6{ + M i au r a h 2 + w aj r a e{. 

Again, as in Case II, a change of frame 63 = R3393 + 4>3Tt leads to the simpli- 



4.4 The proof of Theorem 14. II 



38 



fication of (|4.77cp to 

del = tc<* A ^ + Ml r n a A 6{ + M l au n a A 9% (4.78) 

+ \c l rs 9\ A 6{ + Cl u 61 A k + CI, 6{ A P 3 + \ci v 9% A k + C l uj 6% A %. 

Finally, just as in the reduction from (|4.72p to (|4.75|) . we use a change of frame 
#3 -> 03 + R13O1 + R23O2 to transform (|4~78]) to 

d6\ = \Bi^^ (4.79) 

+ ^ A 6{ + Cl u 61 A % + £ A 0| + ^ A 91 + C l uj 6% A %. 

Equations (|4.77ap , (|4.77b|) and (|4TT9"1) give the desired result. We note once 
more that the structure functions C*„ in (|4.77c[) and (|4.78p are not necessarily 
constant but they are constant in (|4.79[> . 

The reduction of the structure equations for a general p-step solvable algebra 
follows this construction and can be formally defined by induction on p. 

Case IV: Finally, we consider the case where Oess is a generic Lie algebra. 
We use the fact that every real Lie algebra g admits a Levi decomposition into 
a semi-direct sum g — s © r, where r is the radical of g and s is a semi-simple 
subalgebra of g. The radical r is the unique maximal solvable ideal in g - the 
semi-simple component s in the Levi decomposition is not unique. The dual 
space g* then decomposes according to 

0* = ann(r) ann(s). (4.80) 

The fact that t is an ideal implies that 

dann(r) C A 2 (ann(t)) and rfann(s) C g* (g)ann(s). (4-81) 

If M : g — > g is a derivation, then M preserves r and hence M* : ann(r) — * ann(r) . 

Now choose, by a constant coefficient change of basis, 1-forms {61,62} 
adapted to the decomposition (|4.80p with ann(t) = span{ 6\ } and ann(s) = 
span{ 6\ }. In view of (|4.8ip . the structure equations (|4.50p become 

d6[ = ±A r ab 7T a A 7r h + 1^ 7f a A ^ + l -Cl t 9{h6\ + M r as K a A 0> x , (4.82a) 

dOl = \AU *° A + \Bi TT a A 7r' 3 + \c) k 6\ A 9 h 2 + & rj 9\ A 6\ (4.82b) 
+ M i ar * a A6\ + MiA a h6{. 
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In these equations the indices r, s, t range from 1 to sq — dim(ann(r)) and i, j, 
k range from sq + 1 to s — dim(g). 

Since the structure constants CJ t are those for the semi-simple Lie algebra s, 
we can return to the analysis presented in Case I to prove that there is a change 
of coframe §1 — 9{ + <t> r a : K a (see (I4.57[l ) which reduces the structure equations 
(j4~82)) to 

d9[ = -B^p 7r Q A tt^ -+- ~C r st 6{ A 6{ , and (4.83a) 

dO l 2 = \A\ b rt a A n b + tt« A ^ + \c) k 6{ A 9 k + & rj 9{ A 9{ (4.83b) 

+ M: ir n a A9 r 1+ AP a ^ a A9i. 

The structure functions CJ t , C^ k and in (|4.83p are identical to the corre- 
sponding structure constants in (I4.82p . One now checks that the 9\ terms in 
(|4.83b[) do not effect the arguments made in Cases II and III. Thus, by a change 
of frame §2 = R62 + <P^, one can reduce J4.83[) to 

d6{ = -B r af3 w a Aft? + -Cl t 6iA§l, and (4.84a) 

dk = \B l a p * a A^ + ~q k 91 A 9 k + & rj 6\ A 9{ + AP ar ^ A 9\. (4.84b) 

Finally, by a now familiar computation, one sees that a change of frame 62 — * 
62 + R6\ allows one to eliminate the 7r a A 9\ terms in (|4.84b[) . 

This completes our derivation of the structure equations (|4.48|) and the proof 
of Theorem 14. II is. at last, finished. 

Remark 4.9. We collect together a few additional formulas which will be 
needed in the next section for the proof of Lemma 15.81 and the construction 
of the superposition formula. If the forms 9 l and 9 l in (|4.47p satisfy (|4.48p and 
(|4.49p . then it easy to check that coefficients Rp Rp <p l — (/) l a ^ a and ip l — ^ l a 'k a 
satisfy (see (|4.5ip ) 

dR) = -C\ k R) 4> k - R\ Mj, dR) = C\ k Rfy* - R^Nj, (4.85a) 

R\ B * a p = G l a p, R]E £ ab = H l ab , (4.85b) 

d<P = \c\ k #> A<t> k - X -R\A^ dr = ~q k ^Ai> k -^R\F e . (4.85c) 
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The computations leading to (|4.85[) also show that the matrices and 
define automorphisms of the Vessiot algebra, that is 

R-l Cjk = Qm Rj R-k ) an< i Hi Cjk = Rj ^k • (4.86) 

Finally, a series of straightforward computations, based on (|4.44|) . (I4.45|) . (I4.85ap 
and (|4.86[) . shows that the matrix 

A = RQRT 1 satisfies (4.87) 
X\C l jk = C\ m X l j X k n and dX) = C\ m Xf 8 e + X\ C"V V m (4.88) 

I 

We conclude this section by determining the residual freedom in the de- 
termination of the 1-forms 9 l . Specifically, we compute the group of coframe 
transformations which fix the forms 7r a and 7r" and which transform the 6 l by 

¥ = A}0* + a\ where a 1 = S l a Tr a (4.89) 

in such a manner as to preserve the form of the structure equations (|4.48j) . 

If we take the exterior derivative of (|4.89j) and substitute into the structure 
equations for B l we find first, from the 7r a A Q % and the 7r a A 7r a terms, that 
A*-,S* € Int(F) and then that 

dk) = Cl k k]a\ A\C] k = C\ m A\K^, do^\c) k a^a\ (4.90) 

To integrate these equations, let G be a local Lie group whose Lie algebra is 
the Vessiot algebra Dess with structure constants Cj fe . On G construct a coframe 
rf of invariant 1-forms with structure equations dr\ % = \C^ k Tf /\i] k . Then there 

A 

exists ([13 1 Proposition 1.3) a map a: Int(V) — > G such that a 1 = a*{rf). 
Define S(I) = Ad(cr(J). Then 

dS) = G\ k o l S) and d(AK^ (_1) ) ■) = 0. 

and hence the general solution to ()4.90|) is 

a 1 = a*(rj l ), A) = A\ Sj, S = Ad(a(/)), (4.91) 

. 

where the matrix A| is a constant automorphism of the Vessiot algebra. When a 
is a constant map, S is constant inner automorphism so that as far as the general 



solution (|4.9ip is concerned, one can restrict the K\ to a set of representatives 
of the group of outer automorphisms of the Vessiot algebra. 
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Remark 4.10. The bases for the infinitesimal Vessiot transformation groups 

A A V v V . A 

T = { } and r = { Xi } are related by Xj = XIX i, where A is the matrix 
(|4.87p . Let X r , r = l...m, be a basis for the center of T and pick a comple- 
mentary set of vectors X u , u = m 4- 1 . . . s, which complete the X r to a basis 
for T. Choose a similar basis for T. Then, because A defines a Lie algebra 
automorphism, we have that 

v . A v ^ A . A 

V r = A r V t and l u = A U X„ + A u X t , 

where i = 1 ... m and w = m+ 1 . . . a. The second equation in (|4.88[) now implies 
that d\ r s = 0. Consequently, we may always pick bases for the infinitesimal 

A v 

Vessiot transformation groups T and T so that the infinitesimals generators for 
the center coincide. In turn, this implies that any vector in the center of either 
infinitesimal Vessiot transformation group is a infinitesimal symmetry of both 
V and V. I 
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5 The Superposition Formula for Darboux Pairs 

5.1 A preliminary result 

The proof of the superposition formula for a Darboux pair will depend 
upon the following technical result concerning group actions and Maurer-Cartan 
forms. 

Theorem 5.1. Let M be a manifold and G a connected Lie group. Let tu l L and 
u> R be the left and right invariant Maurer-Cartan forms on G, with structure 
equations 

dui = )f i jk u i L Aw k L and dJ R = - l -C) h <J R A u k R . (5.1) 

Let X^ and X^ be the dual basis of left and right invariant vector fields on G. 

[i] Suppose that there are right and left commuting, regular, free group actions 
of G on M, 

p:GxM^M and Ji: G x M -> M, (5.2) 

with common orbits. Denote the infinitesimal generators of these actions by 
(hUXh = X t and faMX?) = Xi. 

[ii] Suppose there are 1-forms u l and lo 1 on M such that [a] Cj l (Xj) = Sj and 
u) l (Xj) = Sj; [b] span{ib 1 } = span{ lo 1 } with ci> l (xo) — w 4 (x ) at some fixed 
point of Xo G M and; [c] 

d£o % = hj) h & A Co k and du l = Af)k & A to k . (5.3) 

Then about each point xq of M there is an open G bi-invariant neighborhood U 
of M and a mapping p: U — > G such that 

p*(w l L )=u)\ p*{uj l R )=uj\ p(jl(x,g)) = p(x) -g, p{fi{g, x)) = g ■ p(x). 

(5.4) 

If M/G is simply connected, then one can take U = M . 

Remark 5.2. This theorem can be viewed as a simple extension (or refinement) 
of three different, well-known results. First, suppose that [i] holds. Then, 
because the actions p and p are free and regular, we may construct p and p 
invariant sets U and U and maps 



p:U->G and p:U^G 



(5.5) 
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which are G equivariant (see, for example, |12j ) with respect to the actions p, 
and p, respectively. Theorem l5.1l states that these maps can be chosen such that 
they are equal on a common bi-invariant domain, p*(iv l L ) = ai 1 and p*(u R ) = uo l . 
Secondly, assume [ii] . Then a fundamental theorem in Lie theory (see Griffiths 
[IT] or Sternburg [53], page 220) asserts that there always exists maps : M — » 
G and <p : M — > G such that 

0*(u i L )=u* and 0*(wjj)=S < . (5.6) 

From this viewpoint, Theorem 15.11 states that these maps can be taken to be 
equal and equivariant with respect to the actions p, and p. Thirdly, suppose 
that just one of the actions in [i] is given, that G acts regularly on M and the 
corresponding forms (say ui l ) in [ii] are given. Then 7r: M — > M/G may be 
viewed as a principal G bundle with (by [ii]) flat connection 1-forms Cj 1 . Then 
Kobayashi and Nomizu [18 (Theorem 9.1, Volume I, page 92) state that there 
is an open set U in M/G and a diffeomorphism ip : 7r~ 1 (L/) — > U x G such 
that ijj*(u) l L ) — Co % (where the Maurer-Cartan forms ui l L define the canonical flat 
connection on U x G). Moreover when M/G is simply connected, then M is 
globally the trivial G bundle with canonical flat connection. Compare this result 
with Corollary 15.61 to Theorem 15. II 

We establish Theorem 15.11 with the following sequence of lemmas. 

Lemma 5.3. Under the hypothesis of Theorem \5.1\ there is a bi-invariant, con- 
nected neigborhoodU around xq and a mapping p: U — > G such that p*(u> l L ) = Cj 1 
and p*{oJ\j) = u) 1 ■ The map p is uniquely defined for a given domain hi. 

Proof. As is customary, we suppose that io l L {e) = io l R {e), where e is the identity 
of G. In accordance with the aforementioned Theorem 9.1 in [18] there is a fx 
invariant open set U in M and a unique smooth map p: U — > G such that 

p*{uj l L )—uj l and p(xq) = e. (5-7) 

Since the p and p orbits coincide, the set ti is actually bi-invariant. Because 
G is connected, we may assume that hi is connected. The uniqueness of p is 
established in Theorem 2.3 (page 220) of [23] , 

In order to prove that p satisfies p*(ui R ) = ui 1 , let A*-(<7) be the matrix for 
the linear transformation 



M*{g)=RloL^ 1 :T:G^T:G 
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with respect to the basis Cj l {e). Then an easy computation shows that 

c4(s)=Aj(<7)<4(<?) and A*( 5lff2 ) = A}( 52 )A^( 5l ). (5.8) 

The Jacobian of the multiplication map m: G x G — > G may be computed in 
terms of A* as 

m*(wi)(5ifl2) = A}( 92 )4( 51 ) +a > i07a). (5.9) 

We take the Lie derivative of the first equation in (|5.8p with respect to Xjf to 
find that L x lK 1 ^ = G£ £ A^ and hence 

dA} = C£,Aju;|. (5.10) 

Similarly, by [ii] [b] , we may write Cj % — AJu^ and calculate from (|5.3p 

dA} = C^Ajw fc . (5.11) 

The combination of equations (|5.7|) . (15.10P and (|5.11|) (and the fact that A* (e) = 
Aj(a;o)) now shows that 

d(p*(A)){\- 1 y k ) = and hence p*(A}) = A*. (5.12) 

(Here we use the connectivity of IX). The equations p*(ui l R ) = lq l now follows 
immediately from (pTT|) . (|5?8|) . and (|5TT2]) . | 

We note, for future use, that (|5.1ip also implies 

L k p = C^{\-% u k + (X-'iC&tf = 0. (5.13) 

This shows, by virtue of the connectivity of G, that 

(5.14) 

Lemma 5.4. With p :IA — > G as m Lemma \5.3\, the maps (p = po p, xa : G — > G 
and ft = po p xo : G — > G satisfy 

£*K)=<4 and $*{lj r )=uj r (5.15) 

and hence (p = tp = idc- 

Proof. The definitions [p x )^(Xl) = Xi and (Ji x )*{Xf) — Xi, together with 
[ii](a), imply that p* Xo {ui l ) — cu l L and p,* (u> 1 ) = w R . Equations (|5.15p then 
follow directly from Lemma l5.3l These equations imply that (p and (/5 are trans- 
lations in G and therefore, since G is connected and <p(e) = <p(e) = e, we have 
that (p = ip = idc- I 
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Lemma 5.5. The maps p: hi — > G satisfies p{p{x, g)) = p(x)-g and p(p,(x, g)) = 
g ■ p(x) for all x £ U and g G G. 

Proof. To prove these equivariance properties of p, define maps 

p = Rg-i o p o p g : M — ► G and p = L„-\ o po p g : M —> G, 

where p g (x) = p(g, x) and p g (x) = p(g, x). We wish to prove that p = p = p. 
Lemma shows that p(xo) = p{xo) =eso that it suffices to show that 

p*(u i R )=u> i and p*(u i L )=u i (5.16) 

and then to appeal to the uniqueness statement in Lemma 15.31 The first equa- 
tion (|5.16|) follows from the fact that the forms lo r are right invariant, the 
second equation in (|5.4[) and (|5.14[) . The proof of the second equation in (|5.16[) 
is similar. | 

Corollary 5.6 (Corollary to Theorem 15. ip . Let S be the submanifold of M 
defined by S = p^ 1 (e), where p is the map (|5.4|) . Then the map 

^■.U^SxG defined by V(x) — (fi(p(x)~ 1 ) x), p(x)) 

is a G bi-equivariant diffeomorphism to an open subset of S x G and satisfies 

*o M*i)=Xi> **{u R )=W- (5-17) 

If M/G is simply-connected, then ^ : M — ► S x G is a global diffeomorphism 

Proof. It is easy to check that sends the intersection of U with the domain 
of p, in M into S x G. To check the equivariance of it is convenient to write 
p{g,x) = x * g and p(g,x) = g~* x. Because the orbits for the two actions 
coincide, we know that for each g S G there is a g' such that g~* x = x * g'. 
We then use (15.41) to calculate 



ty(x% g) — [{x^* g)% p(x% g)^ -1 ', p{x% g)) — {x\ p(x)^ 1 \ p(x)% g) 
= ^l{x) * g, and 

^(g^* x) — ((.9"* x)% p{g~* x)^ 1 ^, p(g~* x)) — (x% g')^* p{x% g')^ 1 ^, g~* p( 

— (fx** p(xp^ 1 ',g"* p{x)) =g*ty(x), 
as required. | 
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5.2 The Superposition Formula for Darboux Pairs 

We are now ready to establish, using the 5-adapted coframe of Section 4.4 
and the technical results of Section 5.1, the superposition formula for any Dar- 
boux pair. 

By way of summary, let us recall that our starting point is the local (0- 
adapted) coframe {9, fj, <x, fj, a} which is adapted to a given Darboux pair 
{ V, V } in the sense that 

/ = span{0, fj, fj}, V = span{0, <x, fj, fj}, V — span{0, fj, <r, fj}. (5.18) 

In our first coframe adaption we adjusted the forms fj and fj to the form 

fj = dJ 2 + R<x and f) = di 2 + Rv. (5.19) 

The fourth coframe adaptations {Ox, 0", fj, <r, f]} and {9y, a, fj, cr, fj} al- 

A v 

lowed us to associate to the Darboux pair { V, V } an abstract Lie algebra 
tiess(V, V). We have 9 l x = QjO Y and we can write 

V(1V = span{0 x , fj, fj}, (5.20) 

From the 4-adapted coframes we then constructed the final coframes 

8 l = Rfy + ^ a , 6 l = R^y + ^ a 7r a , (5.21) 

where ft = [a, fj], and tt — [er, fj], with structure equations 

do 1 = l -G^^+ l -c) k e^e k , d¥ = att 6 -\c] k v A§ k . (5.22) 

The coefficients G^a are functions of the first integrals i while the H l ab are 
functions of the / . The two 5-adapted coframes 

{9, <r, fj, <j , fj} and {9, a, fj, a, fj} (5.23) 

A _____ 

are not adapted to the Darboux pair { V, V } (in the sense of (|2.6p ) although 
we do have 

V = span{ 9, <r, fj, fj} and V = span{ 9, fj, cr, fj}. (5-24) 

Definition 5.7. Let { V , V, } be a Darboux pair with 5-adapted coframes (|5.23[) . 
Let { Xi, V a , V a } and { Xi, W a , W a } be the dual frames corresponding to the 
5 adapted coframes so that, in particular, 

e\x 3 ) = 5), ^(Xj) = o, e^Xj) = 5), n a {x 3 ) = o. (5.25) 
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The infinitesimal Vessiot group actions for the Darboux pair { V, V } are 
defined by the Lie algebras of vector fields Xi and Xi. These define (local) 
groups actions 

£L:GxM^M and jl: G x M M. (5.26) 
We take ft to be a right action on M and jl to be a left action. 

The vectors fields X { and X { are related by (see (|4.87j) and (|4.88|) ) 

Xi = XiXj and satisfy [Xi, Xj } = 0. (5.27) 
From the structure equations (|5.22p we find that 

C x P = &J k and £x/ = q k 8 k . (5.28) 

= £x 3 * a = £xfi a = ^ 3 7f Q = 0. (5.29) 

The actions /t and fi commute because [Xi,Xj] = 0. In what follows we shall 
assume that these actions are regular. 

Now let G be a Lie group with Lie algebra vzss(V ,V). Let us % L and be 
the left and right invariant Maurer-Cartan forms on G, with structure equations 
(|5.ip . We shall assume that these two coframes on G coincide at the identity e. 
Let X^ and Xf- be the dual basis of left and right invariant vector fields. 

To apply Theorem 15. II we need to identify the forms u? and ui i . 

Lemma 5.8. The forms Cj 1 and oj l defined in terms of the 5-adapted coframes 
($MB or (|5^TD ) by 

Cj* = + Xiipirr" and d* = ¥ + a)^ a Tt a , (5.30) 

where X = RQR and fi = A (see Remark have the same span and 
satisfy the structure equations 

du l = hj) k Cj j A Co k and du l = Af)k & A uj k . (5.31) 

Proof. To check that the forms Cj 1 and uj 1 have the same span we simply use the 
various definitions given here to calculate 

Co 1 = Q l + X) ipi 7f Q = R) 4 + <P a Z a + X)^ a 

= r) q{ e k Y + tfjr + \)^ a = (rq it 1 )) Ri e k Y + ^J a + x) ^ tP 

= X) (Ri 9 k Y + ^ + (X-% 0* n a ) = A* (8 j + <f> k n") = \*&. (5.32) 
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We now calculate dCj 1 using flOSJ), P~g5]) . P~g6")) and (|Q5|) to find that 
d£j l = d§ 1 + dX) A V j + A} A 

- ^ *° A ^ + ~Cj fe # A 9 k + (Cj m X™ £ + V> m ) A ^ 

+ A} (-^ fe ^ A ^ fe - ^ 7f° A *") 

= l -C) k frfx6 k + C\ m A™ £ A ^ + |A| C£ fe ^ A ^ = iq fc ^ A Q k . 
The derivation of the structure equations for uj 1 is similar. | 

Let xq € M. By a change of frame at xq we may suppose that lo 1 (xq) = 
ui 1 (xq). At this point all the hypothesis for Theorem 15.11 are satisfied and, in 
accordance with this theorem and Corollary [SHI we can construct maps 

p:U^G and f:M^5xG (5.33) 

which satisfy (|5.4p . (|5.17|) and p(xo) = e. 

Define S\ and 52 to be the integral manifolds through xo for the restriction 
of V (Qc) and y (oo) to S. Because 

T*S X = t>i oo) © V^* for all x E S, (5.34) 

we may choose 5i and £2 small enough so as to be assured of the existence of 
a local diffeomorphism 

X :5 1 x5 2 ^5 , (5.35) 
where Sq is an open set in S. The map \ satisfies (with s = x( s 1t s 2) and 

x * ( t>(°°))=T s ; 4 S 2 and x*^) = T* S x . (5.36) 

A v A v 

Set xO^i) ^2) = x o- Hence, for any first integrals I and I of V and y, we have 
I(x(si,s 2 )=I(x(x° 1 ,s 2 )) and 7( X (si, s 2 )) = I( X (s u x%)). (5.37) 
Put W = $ _1 (<So)- 

To define the superposition formula for the Darboux pair {V, V}, we take 
Mi to be the (maximal, connected) integral manifold of y*- 00 - 1 through xo in Uq 



5.2 The Superposition Formula for Darboux Pairs 



49 



and M 2 to be the (maximal, connected) integral manifold of y K °°' through xq 
in IAq. The group actions p and p and the maps p and ^ all restrict to maps 

Uj : G x Mi -> Mj, fa: G x Mi -> Mi, 
^ ' M (5.38) 

Pi : Mi -> G, and : Mi ^5,xG 

for i= 1,2. 

With the maps 4" and x defined by (|5.4[) and (|5.35[) and with these choices 
for Mi and M 2 , we define 

E : Mi x M 2 — > M by Ufa, x 2 ) = m(<?i • 92, s 2 )), (5.39) 

where € Mi, and \&(:Ei) = (si,gi). 

To prove that E defines a superposition formula, we must calculate the 
pullback by S of an adapted coframe on M. To this end, we introduce the 
following notation. For any differential form a on M, denote the restriction of 
a to Mj by a^. If / is a function on M and Xi S Mj, then /(&») denotes the 
value of f at Xi, viewed as a point of M . We also remark that the equivariance 
of p with respect to the group action p and p implies that E satisfies 

(p o E)(xi,x 2 ) = p(p(gi ■ 92, x(si, s 2 )) = gi ■ 92 ■ p(x(si, S2)) 

(5.40) 

= 9i-92 = Pi(xi) ■ p 2 (x 2 ) = (mo (pi x p 2 ))(a;i,x 2 ). 

Lemma 5.9. TTie pullback by E of the ^-adapted coframe (see Theorm \4.6\ ) on 
M is given in terms of the restrictions of the 5-adapted coframes (|5.23|) on Mi 
and M 2 by 

E*(«x) = <r 2 , E*(<x)=*l> X*(f))=V2, E*ft) = tfe, (5.41) 
E*(i?6» x ) = A(0 x + 2 ). (5.42) 

Proof. We begin with the observation that the first integrals / and / are 
invariants of p and therefore, by (|5.37[) . 

E*(I a )(xi,x 2 ) = I a (p(gig 2 ,x(si,s 2 ))) = I a (x(si,s 2 )) 
= I a (x(xoi,s 2 )) = I a {p{g 2 ,x{x i,s 2 ))) 

= / (x 2 ), and similarly (5.43) 

Z*(I a )( Xl ,x 2 ) = I a ( Xl ). (5.44) 
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For future reference, we note (see (|5.2ip ~) that 

Y?{G%){x l ,x 2 ) = G%{x 1 ) and x 2 ) = Gfofa). (5.45) 

From equations (|5.43|) and (|5.44|) it follows directly that 

£*(<t) = <t 2 , E*(*)=*i, (5.46) 
T,*(rj) = E*(dl + R<x) = rj 2 and S* (fj) = S* (dl + R <x) = ^ . (5.47) 

Equations (|5.4ip are therefore established. 
Equations (|5~2Tj) and (fOP]) imply that 

^4=^-^7r°-AiV^ a . (5.48) 

A. A 

We use this equation to calculate H*{Ra6 3 x ). 

Since € Int(V') and i// a S Int(V"), equations (f5T43|) and (f5T44|) imply that 

E*(0 4 Q ^)( a ; 1 ,a ;2 ) = 4 Q (x 2 )^ and E*(V£ Tr")^,^) = ^(n) tt?. (5.49) 

We deduce from equations (fST8j) . ([5321) and ijOOf that 

£*(Afc)(a;i,ar2) = £*(p*A l fe )(xi, ar 2 ) = A l k ( gi ■ g 2 ) = 

(5.50) 

= A}(g 2 )Ai( gi ) = X}(x 2 )4( Xl ). 

To calculate £*(a) J ) it is helpful to introduce the projection maps tt, : G x 
G — ► G and to note that 7Tj(pi x p 2 ){xx, x 2 ) = g%, i = 1, 2. We can then re-write 
([S3)) as 

m*K)=^(AjKK)+^K). 
This equation, together with (|5.4|) and (|5.40p leads to 

S*((i l )(a;i,a;2) = S*(p*(w^ ac 2 ) = (p* x paXm*^))^, z 2 ) 

= (Pi x P 2 *)(7rI(A}) 7 r 1 *K) +7r 2 *(wi))( a!l , a : 2 ) = ( 5 ' 51 ) 
= A}(x 2 )w^ +tj 2 . 

Finally, the combination of equations (I5.32[) and (|5.48p - (|5.51|) allows us to 
calculate 

£* (^ P x ) = £* - # n a - XI ^ tt<*) 

= X)(x 2 ) Cj{ +u>1- J>i(x 2 ) ^ - X)(x 2 )Xi(xi)ip k M *? 

= X){x 2 ) (# - Xi (zOV^i) X? ) + A}(z 2 ) (4 - n{{x 2 ) <f> k a (x 2 ) ^ ) 

= X)(x 2 )ei + X)(x 2 )ei, (5.52) 
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which proves (|5.42[) . | 

Recall that for any Darboux pair { V, V }, the EDS V ffl V is given by Defi- 
nition 12.51 

A V 

Theorem 5.10 (The Superposition Formula). Let {V, V} define a Dar- 
boux pair on M. Define Mi and M 2 as above and let W\ and Wi be the restric- 

— /\ 

Hons ofV and V to Mi and M2, respectively. Then the map S: Mi x M2 — > M, 
defined by (|5.39|) . satisfies 

Z*{VnV) C Wi + W 2 and E* (V ffl V) C Wi + VV 2 . (5.53) 

Thus E defines a superposition formula for V PR V with respect to the Pfaffian 
systems Wi and W 2 ■ 

Proof. To prove this theorem we need only explicitly list the generators for 
Wi + W 2l Wi + W 2 , VC\V and V ffl V and check that, using Lemma |[5TT7|) . that 
the latter pullback into the former by E. For this we use the definitions of V 
and V in terms of the 4-adapted and 5-adapted coframes bv (|5.20p and (|5.24[) . 

From the definition of Mi and M 2 as integral manifolds of V^ 00 ' and , 
we have that 

fj 1 = 0, <T! = on Mi and fj 2 = 0, <r 2 = on M 2 

and hence the two 5-adapted coframes on M restrict to coframes 

{ §1, fji, & 1 } for Mi and { 2 , fj 2 , ^2 } for M 2 . (5.54) 

These naturally combine to give a coframe on Mi x M 2 . The 1-form generators 
for the Pfaffian systems Wi and W 2 are then given by 

VFi=span{#i, } and VF 2 = span{ 2 , j} 2 }. (5.55) 

The differential system Wi + W 2 is algebraically generated by the 1-forms 
(|5.55|) and their exterior derivatives (see (|2. 13|) . (I5.22| ). that is, 

Wi + W 2 = {9i, i)i, 6 2 , fj 2 , dfji, dfj 2 , dOi, d6 2 }. 

The restrictions of (15.221) to Mi and M 2 lead to the structure equations 

d0\ = -Gi7Ti A 7Ti mod {©i} = -G1CT1 A <ti mod {#1,771} and 

d0 2 = \ i H 2 -k 2 A 7r 2 mod {6> 2 } = \-H 2 & 2 A <r 2 mod{0 2 ,f) 2 } 
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from which it follows that Wi + W2 is algebraically generated by 
Wi+W 2 = {01, »7i, 02, i? 2 , FrfiA&x, F 2 &2 A & 2 ,H 2 &2 A o- 2 }. (5.56) 

The Pfaffian system FnT^ is generated by the forms {0, fj } 17} or cquivalcntly 
(see (|5.20[) ) by {Ox, V, v}- The first inclusion in (|5.53[) follows immediately from 
the Lemma [531 and (|5.55|) . 

A, 

By definition, the differential system V ffl V is algebraically generated, in 
terms of the 1-adapted coframe (see Theorem I2.9|) . as 

VfflV = {0, r), rj, di), di), A&A&, B a A <r } 

or, equally as well, in terms of the 4-adapted coframe by 

VfflV = {0, r), 17, F&Aa, F&A&,A& A <t, BaAa}. 

In this latter equation the coefficients A and B are those appearing in (|4 . 25[) . 
By (|4.46|) and (|4.85[) . we also have that 

V ffiV = {R8 X , r), 17, F&A&, Fa Aa,H & A&, G&Aa}. (5.57) 

With the algebraic generators for V ffl V in this form, it is easy to calculate the 
pullback of V ffl V to M x x M 2 by E. By using (|5~4Tj) . (|5~l2")) . (|532")l we find that 

E*(V ffl V) = {0i+02, 172, *h, Ft&xA&i, F20-2A6-2, if 2 <5- 2 A<T2, GiffiAffi} 

(5.58) 

and the second inclusion in (|5.53[) is firmly established. | 

One may use the first integrals of V and V and the map p: U — ► G, and 
local coordinates on G to define local coordinates 

(7 a = I a {x), 7 a = z l = p\x)) (5.59) 

on M and induced coodinates (J",zJ) on Mx and (-f^i^I) 011 -^2- In these 
coordinates the superposition formula becomes 

E((/?,4), (3,4)) = (/" = ^ = = (*x) ■ (4))- (5.60) 

We shall use this formula extensively in the examples in Section 6. 

To prove that the superposition map E is subjective, at the level of integral 
manifolds, we use the concept of an integrable extension of a differential system 
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X on M [5]. This is a differential system J" on a manifold N together with a 
submersion <p: N — > M such that ip*(X) C c7 and such that the quotient dif- 
ferential system J /f*X is a completely integrable Pfafhan system. This means 
that there are 1-forms d l <E J such that J is algebraically generated by the 
forms { -& e } U <^*T and 

d^ = mod { 0* } U p*I (5.61) 

or, equivalently, modulo the 1-forms in J and the 2-forms in ip*T. Under these 
conditions the map tp is guaranteed to be a local surjection from the integral 
manifolds of J to the integral elements of X. Indeed, as argued in [5], let P C M 
be an integral manifold of X defined in a neighborhood of x E M. Choose a 
point y £ N such that (p(y) = x. Then P = tf~ 1 (P) is a submanifold on N 
containing y and the restriction of J to P is a completely integrable Pfaffian 
system J. By the Frobenius theorem (applied to J as a Pfaffian system on P) 
there is locally a unique integral manifold Q for J" through The manifold Q 
is then an integral manifold of J which projects by ip to the original integral 
manifold P on some open neighborhood of x. 

Corollary 5.11. Let {V, V} be a Darboux pair on M. Then the EDS Wi +W 2 

A V 

on Mi x M2 is an integrable extension of V ffi V on M with respect to the 
superposition formula S: Mi x M2 — > M. 

Proof. In view of (|5 . 56[) and (|5 . 58[) we may take the differential ideal generated 
by T — { 6{\ as a complement to S*(D ffl V) in Wi + W 2 - The 1-forms 0i are 
closed modulo 6\, rj 1 and G1&1 A <ti and therefore the generators for T are 
closed, modulo T and modulo the 1-forms and 2-forms in E*(V Ffi V). Hence 
Wi + W2 is an integral extension of V ffi V. | 

Corollary 5.12. LetX be a decomposable, Darboux integrable Pfaffian system. 
Then the map S : Mi x M2 — ► M is a superposition formula for X with respect 
to the EDS Wi + W2 which is locally surjective on integral manifolds. 

Proof. If { V, V } is the Darboux pair for X, then X = V FR V and the corollary 
follows from Theorem 15.101 and Corollary 15.111 | 
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5.3 Superposition Formulas and Symmetry Reduction of 
Differential Systems 

In this section we shall use the superposition formula established in Theo- 
rem 15.101 to prove that if X is any decomposable differential system for which 

A v 

the singular Pfafhan systems {V, V} form a Darboux pair, then X can be real- 
ized as a quotient differential system, with respect to its Vessiot group, by the 
construction given in Corollary [331 

To precisely formulate this result, it is useful to first summarize the essential 
results of Sections 4 and 5.1-5.2. We have constructed, through the coframe 
adaptation of Section 4, a local Lie group G and local right and left group actions 
p, p, : G x M — > M, For the sake of simplicity, let us suppose that G is a Lie 
group and that these actions are globally defined. The infinitesimal generators 
for these actions are the vector fields Xi and Xi, defined by the duals of the 
5-adapted coframes. 

[i] The actions p and p are free actions with the same orbits (the vector fields 
Xi and Xi are pointwise independent and related by Xi — XjXj). 

[ii] The actions p and p commute (see (15.271) ). 

[iii] The actions /i and p are symmetry groups of V and V , respectively (see 

[iv] Each integral manifold of V°° or V"°° is fixed by both actions p and Jl (see 

[v] We have defined i\ : M\ — > M and ii : Mi — > M to be fixed integral 
manifolds of V and V and W\ and W2 to be the restrictions of V and 
V to these integral manifolds. 

[vi] Properties [i] [iv] imply that the actions p, and Jl restrict to actions on pi 
and pi on Mj. These restricted actions are free, pi and /i2 are symmetries 
of W\ and W2 respectively. The actions pi and /i2 are transverse to Wi 
and W 2 (see (IQ51 V 

The diagonal action 8: G x (Mi x M 2 ) — > Mi x M 2 is defined as the left action 



S(h,(xi,x 2 )) = (pi(h 1 ,xi), p 2 (h,x 2 ))- 



(5.62) 
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Note that the infinitesimal generators for 5 are 

Z i = -X li +X„. (5.63) 

Granted that the action 8 is regular, we then have that all the hypothesis of 
Corollary 13.41 are satisfied. We can therefore construct the quotient manifold 
q : Mi x M 2 -> (Mi x M 2 )/G, the quotient differential system J = (Wi +W 2 )/G 
and the Darboux pairs { U, U} (see p. 171) ). 

We use the superposition formula (I5.39[) to identify the manifold M with 
(Mi x M 2 )/G and the original differential system I with the quotient system 
J. 

Theorem 5.13. LetX be a decomposable differential system on M whose sin- 
gular Pfaffian systems { V, V } define a Darboux pair. Let G be the Vessiot 
group for { V, V } and define Pfaffian systems W\ on Mi and W 2 on M 2 as 
above. Then the manifold M can be identified as the quotient of Mi x M 2 by 
the diagonal action S of the Vessiot group G, the superposition formula E is the 
quotient map, and X — (Wi + W 2 )/G. 

All the manifolds, actions, and differential systems appearing in Theorem 
15.131 are presented in the following diagram: 

Proof of Theorem \5.13\ The following elementary facts are needed. 

[i] The superposition map E : Mi x M 2 — » M is invariant with respect to the 
diagonal action S of the Vessiot group G on Mi x M 2 . 

[ii] For each point x — (xi,x 2 ) £ Mi x M 2 , kerS*(a;) = kerq^x). 

[iii] E*(V) = [Wi 8 A 1 (M 2 )] sb and E*(V) = [A 1 (Mi) 8 W 2 ] sh . 

Facts [i] and [ii] show that E : Mi x M 2 — > M can be identified with q : Mi x 
M2 — > (Mi x M 2 )/G. Fact [iii] shows that V =U and V =14 (in the notation 
of Corollary[3!]) so that 1 = V ffl V = W ffl W = J . 

To prove [i], let x\ £ Mi, X2 € M2 and h £ G. If ^(xi) = (si,ffi) and 
^(^2) = (s 2 ,g 2 ) then it is a simple matter to check, using the G bi-equivariance 
of p, that jli(h~ l , £1) = (si, (71 • /i" 1 ) and /X2(/i ; ^2) = (s 2 ,h- g 2 ). The invariance 
of E then follows immediately from its definition (|5.39[) . 

Lemma 15.91 and (|5.63[) show that 

ker E* = span-j^ - dg 2 } — span{-Aij + X 2i \ = span{ Z t } (5.64) 
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which proves [ii] . 

Lemma 15.91 and equations (|5.18[) show that 

E*(l>) - £*({0, «x, fj, 77} = { 0i + §2, & a , V2, fit } (5-65) 

while ([5~5"4]) and ([535]) give 

Wi © A\M 2 ) = { 0!, 77^ 2 , 77 2! 6-2 } (5.66) 

in which case [iii] follows immediately from (|5.64j) . 

To make precise the identification of M with the quotient manifold M = 
(Mi x M2)/G, define a map T : M — > M as follows. For each point x € M, pick 
a point (xi,x 2 ) G Mi x M2 such that q(xi,:E2) = x and let T(x) — £(2:1, Xa)- 
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The diagram 



Mi 



commutes and, on the domain of any (local) cross-section £ of q, one has T = 
So(. This observation and facts [i] and [ii] then suffice to show that T is a 
well-defined, smooth diffeomorphism. Moreover, the same computations used 
to establish fact [i] show that S, and hence T, is G equivariant with respect to 
the action Ji\ on M\ x M 2 and jl on M and also G equivariant with respect to 
the action /t2 on Mi x Mi and jl on M. 

Finally, we recall that U and U may be calculated from the cross-section £ 
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and the 6 semi-basis forms by 

U = C([W 1 ®T*(M 2 )} sh ) and U = C*(T*(Mi) W 2 ] sb ) (5.67) 
in which case [iii] implies that T*(V) = U and T*(V) = U. | 

Remark 5.14. Finally we remark that all the results of this section remain 
valid in so long as the singular systems for J have algebraic generators 

V = {0, 17, 17, <r, P&A&} and V = { 0, 17, f), &, Q&Aa,}. (5.68) 

v A 

where P € Int(V) and Q S Int(V). Such systems need not be Pfaffian. 
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6 Examples 

In this section we illustrate our general theory with a variety of examples. 
Examples 16.11 and 16.21 are taken from the classical literature and are simple 
enough that most of the computations can be explicitly given. We consider PDE 
where the unknown functions take values in a group or in a non-commutative 
algebra in Examples 16.31 and 16.41 In Example 16.51 we present some novel exam- 
ples of Darboux integrable systems constructed by the coupling of a nonlinear, 
Darboux integrable scalar equation to a linear or Moutard-type equation. A 
Toda lattice system and a wave map system are explicitly integrated in Exam- 
ples 16.61 and 16.71 In Example 16.81 we solve some non-linear, over-determined 
systems in 3 independent variables. 

Example 6.1. As our first example we shall find the closed- form, general so- 
lution to 



This example is taken from Goursat's well-known classification (Equation VI) 
of Darboux integrable equations [16J and is simple enough that all the steps 
leading to the solution can be explicitly given. See also Vessiot 28J (pages 
9-22) or Stomark [23] (pages 350-356). 

The canonical Pfaffian system for (|6.ip is / = { a 1 , a 2 , a 3 }, where 

a 1 — du — pdx — q dy, or = dp — r dx — vpq dy, a 3 — dq — vpq dx — tdy 

and v — l/(u — x). The associated singular Pfaffian systems are 

V = {a l , dx, dr — q(vr + v p) dy } and V — {a l , dy, dt — tpvdy}. (6.2) 

A v 

The first integrals for V and V are 

fl t 2 P , , ,2 rl t 2 t la Q\ 

I = x, I = vp, I = vr + v p, 1 = y, I = y (6.3) 

q 

and we easily calculate that V^ 00 ^ n V — { fj 1 } and V n V^ 00 -* = { }, where 

rf =dP + ((f) 2 - P) dl 1 = va 2 - v 2 pa 1 . (6.4) 
A 1-adaptcd coframc (|2.6| is therefore given by 
9 1 = a\ 6 2 = a 3 , a 1 = df , a 2 = dl 3 , fj\ ^ = df , a 2 = dP . (6.5) 
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We relabel the coframes elements by tt 1 — a 1 , tt 2 — a 2 , tt 3 = r) 1 , tt 1 = a 1 , 
tt 2 = <j and calculate 

dn 3 = -2I 2 n 1 AP+n 1 An 2 , 

dO 1 = {u - x) it 1 A n 3 + I 2 n 1 A 6 1 + n 1 A 6 2 , (6.6) 
dO 2 = q n 1 A n 3 + q tt 1 A tt 2 + f tt 1 A 6 2 + P n 1 A 6 2 . 

This coframe satisfies the structure equations (|4.13[) and is therefore 2-adapted. 

The next step, described in Section 4.2, is to eliminate the n a A 6 l terms 
from (|6.6[) . We calculate the distributions U and U (see (|4.15p ) and their derived 
flags to be 

U = {dp + ((I 2 ) 2 P)dp, dp, dp }, U = {dp, dp}, (6.7) 
U {oo) = UU{(u-x)d i +qd g2 , dp}, U (ca) = UU{qd 6 2, qd Ql }. 

and then, in accordance with (|4.f 9[) - (|4.21[) . define 

Xx = (u-x)dgi +qdg2, X 2 = dgi, Y 1 =qdg 2 , Y 2 = dp, (6.8) 

The coframes dual to the vector fields { X,, U, U } and { Yi, U, U } are the 3- 
adapted coframes 

e x = -e 2 , el^o^^^e 2 , e Y = -e 2 , e 2 Y = -e l with (6.9) 
q q q q 

dO x = n 1 A n 3 + tt 1 A n 2 = d6 Y , 

d6 x = -{u - x) tt 1 A tt 2 + 6 X A 9% + n 1 A 9 X + I 2 tt 1 A 9 X , ( 6 10 ) 

d e 2 Y = V—^. n 1 a % 3 - e Y a e 2 + tt 1 a e Y - T n 1 a e 2 Y . 
q 

The coframe (|6.9|) is in fact 4-adapted and hence the Vessiot algebra for (|6.1|) is 
a 2 dimensional non-abelian Lie algebra. 

We may skip the adaptations given in Section 4.3 and move on to the final 
adaptations given in Section 4.4. The Vessiot algebra is 1-step solvable and 
the structure equations (|6.10p are precisely of the form (|4.68[) . The change of 
coframe 1 = 6 X + Pn 1 transforms the structure equations (|6.10[) to the form 
(|4.70[) . The change of coframe 9 2 = 9 X — xO 1 leads to the 5-adapted coframe 
{9 1 ,6 2 } with structure equations 

dO 1 =Tt 1 An 2 , dd 2 = -un 1 Att 2 + 6 1 A6 2 . (6.11) 
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Similarly, the 5-adapted coframe 9 1 = 9\ + IPn 1 , 9 2 = 9\ — yO 1 satisfies 

d$ l d9 2 = U ~ X ~ yq 7t 1 An 3 -0 1 Ae 2 . (6.12) 

q 



Before continuing we remark that the vector fields X\ , X 2 , defined by 

A A 

are given in terms of the dual vector fields X\ and X 2 , computed from the 5- 

< — A A 

adapted coframe, by X\ = X\ — xX 2 and X 2 = X 2 . These vector field systems 
have the same orbits and structure equations but the actions are evidently 
different and it is the latter action that is needed to properly construct the 
superposition formula. 
The forms (jOU)) arc 



9 + I n = — , ui — 9 — ul n = du — qdy — u — , 
Q q 

ui = 9 + 1 it — — , ui = 9 + ( — h yvp) 7r = dy — y — . 

q q q q 



(6.13) 



The Vessiot group for (|6.ip is the matrix group 
forms 



1 b 
a 



with Maurer-Cartan 



, da i b i da 2 db 

ui L = — , u> L =db da, w fl = — , u R = —. 

a a a a 

The map p: M — > G defined by a = q and b = u — yq satisfies (I5~ 

Finally, if we introduce coordinates y\ = y, u\ = u, qi = q, t\ = t on the 
V 1 " 00 ^ integral manifold Mi = { I 1 = 0, I 2 = 0, J 3 = 0} and = u 2 = u, 
p 2 = p, qi = u, r 2 = r on the V integral manifold M 2 = { I = 0,1 =0}, 
then the superposition formula (|5.60p is 

P P2 r p r 2 p 2 



X = x 2 , 



u — x u 2 ~x 2 u — x (u — x) 2 u 2 — x 2 (u 2 — x 2 ) 



2 



y = Vl, - = —, q = qiq 2 , u - yq = u 2 + (ui - 2/iqi)<?2, 

q qi 

or, explicitly in terms of the original coordinates {x, y, u,p, q, r, t} on M, 

u\q 2 

x = x 2 , y = yi, u = u 2 + q 2 ui, p=(H )P2, q = qiq2, 

U2 ~ X2 (6.14) 



r = (H Jf2 + 7 7n, t = tiq 2 . 

U 2 — X2 («2 — ^2) 



It remains to find the integral manifolds for W and W. Restricted to Mi, 

A A 

the Pfaffian system becomes W = { dui — qi dyi, dqi — tidyi } with integral 
manifolds 

yi =0,ui = /(/?), qi = /'(/?), ti = /"(/?). (6.15) 
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Restricted to M2, the Pfaffian system V becomes W = {du2 — P2<ix2, dp2 — 

P2Q2 . v 

r 2(1x2, dq2 dx2}- To find the integral manifolds of W we calculate the 

112 — X2 

^ (2) 

second derived Pfaffian system to be W — { dq 2 du2 } which leads 

U2 — x 2 

to the equation 52 du2 — 1*2 dq2 + X2 dq2 — or 

d (^)_ X2d (l) =0 or d (^l) + l dX2=0 . 

The integral manifolds for W are therefore given by 

x 2 = a, u 2 = x 2 - g(a)/g'(a), q 2 = -l/g'(a), (6.16) 

with P2 and r2 determined algebraically from the vanishing of the first and 
second forms in W. The substitution of (|6.15p and (|6.16|) into the superposition 
formula (|6.14p leads to the closed form general solution 

u= - f{y ]- g{x) +x (6.17) 
9{x) 

for (jBTTj) . 
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Example 6.2. In this example we shall construct the superposition formula for 
the Pfaffian system / = {a 1 , a 2 , a 3 }, where 

a 1 — du — pdx — q dy, a 2 — dp + — r (tan br — br) dx — s dy, 

(6.18) 

a 3 = dq — s dx — b(br + cot br) dy. 

The coordinates for this example are (x,y,u,p,q, s,t) and & is a parameter. 
This example nicely illustrates the various coframe adaptations in Sections 4 
and has a surprising connection with some of Cartan's results in the celebrated 
5 variables paper [TJ. The values b = 1, b = (the Pffaffian system (|6 . 1 8|) 

remains real) and the limiting value 6 = give the three Pfaffian systems for 
the equations 

which are Darboux integrable at the 2-jet level ([2], pages 373-374) and @J, 
pages 400-411). The case b = is treated in Goursat([H] Vol 2, page 130-132. 

One easily calculates V°° n V — V D V 00 = { } and that the first integrals 
for V and V are 



i l =s + Tl I 2 = —{x + 6 2 y)/ 1 + q + b 2 p, 
I 1 = s-t, I 2 = -(x-b 2 y)I 1 +q-b 2 p. 

We immediately arrive at the 2-adapted coframe 

Tt 1 =dl 1 , w 2 = dl 2 , ir^dl 1 , Tr 2 =dl 2 , 9 1 = 2a\ 

6 2 = -bcotbra 2 + -ta,nbra 3 , 9 3 = -boot bra 2 - - tan br a 3 , 
b b 



(6.19) 



(6.20) 



with structure equations 

d6 l = (x + b 2 y) Tt 1 A e 2 + n 2 A 6 2 - (x - b 2 y) tt 1 A 9 3 - it 2 A 9 3 , 

d9 2 = n 1 A n 2 -bcot2brn 1 A 9 2 + bcsc2bTTr 1 A <9 3 , (6.21) 

d6 3 = 7T 1 A tt 2 — b esc 2br n 1 A 6 2 + b cot 2br tt 1 A 9 3 . 

To compute the 3-adapted coframe Q % Xl we simply calculate the derived flag 

A 

for the 2 dimensional distribution U = { d%i, dc^i } (see Section 4.2). From the 
structure equations (|6.21|) we find 

U (1) =UU{ -d 6 2 } and 

A f21 a , (6.22) 

U { > = UU{ -dg2,(x + b 2 y)d g i - 6cot26r9 S 2 - b esc 2br dgs , d g i }. 
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We take { X\, X2, X3 } to be the last 3 vectors in U and calculate 

[X U X 2 ] = b 2 X 3 . (6.23) 

Therefore the Vessiot algebra for the Pfaffian system (|6.18p is abelian if b = 
and nilpotent otherwise. The 3-adapted coframe 6 x (see Theorem 14.4)) is 

°x = -° 2 +cos26r6» 3 , 9 X = -isin25r6» 3 , X = 8 1 + ~(x + b 2 y) sm2b T 8 3 , 
and the structure equations (|4.22[) arc 



dd x = -tt 1 A^ 2 +COS2&T7T 1 Att 2 + b 2 n 1 A6 X , 



d6 x = -isin26r7f 1 A tt 2 - tt 1 A > 



x- 



(6.24) 



1 



d &x = -^(x + b 2 y) sin 267-7^ A n 2 - n 2 A 9 X - b 2 9 x A 9 



This coframe is actually 4-adapted. 

We labeled the vectors in the derived flag (|6.22|) so that the last vector X3 
spans the derived algebra of the Vessiot algebra. By doing so we are assured 
that the structure equations (|6.24|) are of the precise form (|4.68ap - (14.68b|) . The 
analysis at this point refers back to Case II and the structure equations (|4.59|) . 
as applied to just the first two equations in (|6.24|) . The matrices M and R (see 
equations (|4.60p and ()4.6ip ) are found to be 



M 



b 2 * 1 

-TT 1 



and R 



1 1 
cos bl —bsmbl 



sin bl 



cos bl 



(6.25) 



We then compute the 2-forms \ % an d the 1-forms 4> l (see (I4.63p ~) to be 



7 r 1 A 1 . ^ 2 2 

cos bl it A 7r , x 



• sin bl w A tt 



b 2 = --I 2 smbl 1 n 1 , 
b 



(6.26) 



(j) 1 = -I^coabl 1 ^ 1 , 
so that the forms (|4.69|) are 

6\= cos bl 1 6 X -bsinbl 1 9 X - I 2 cos bl 1 , 



I 2 = -sinb^ei + cosbl^l - -I sinbl tt . 
1 b x x b 



(6.27) 



The structure equations (|6.24|) are now reduced to the form (14. 70[) . The final 
required frame change (|4.74[) leads to the 5-adapted coframe 

2 



I, 6 2 = 9 2 , 9^ = 6 x + I 2 cosbI 1 6\ + bI 2 smbre 2 1 + -(I 2 ) 2 K 2 {&.2S) 
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with structure equations 



dO 1 = cos&Ptt 1 Att 2 , d9 2 = - smbl 1 tt 1 A tt 2 , 

o 

d # 3 = - ((at + b 2 y) sin 2br + bl 2 cos 2br) tt 1 A tt 2 - b 2 ^ A 2 . 



(6.29) 



The diffeomorphism $(x,y,u,p,q,s,T) — (x, — y, — g, — s, r) is an invo- 
lution for (|6. 18|) in the sense of Remark 12.131 and can therefore be used to find 
the 5 adapted coframe. 

The forms (jOO)) arc 



Cj 1 = 9 1 +I 2 COsbI 1 fc\ 



1 



r +7 2 -sin67 1 7r 1 , 
o 

1 



cD 3 = 6 3 + An\ (6.30) 



where A = I (I 2 cos2br + ~(x + b 2 y)s\i\2br - -^I 2 )- 

The Vessiot group for this example is the matrix group 



1 z 1 z 3 + b 2 z 1 z 2 /2 
1 6 2 z 2 
1 



(6.31) 



with multiplication 

b 2 

z 1 = z\ + zl, z 2 = zj + z 2 , z z = z\ + z\-—{z\zl-z\z\) (6.32) 
and left invariant forms 



Cj l = dz\ Co 2 = dz 2 , & 3 = dz 3 + ^-(z 1 dz 2 -z 2 dz 1 ) 
The map p (Theorem I5.1|) is then found to be 



(6.33) 



z 1 = — (— 2xsin6s sin&r + 2b 2 y cos bs cos br — bl 2 sin 61 1 + bl 2 sin bl 1 ), 
z = — ( 2fex cos bs sin 6r + 2yb sin 6s cos 6r + I b cos bl —I b cos I ), 



(6.34) 



z 3 = 2u ~ — (sin26r - 26rcos2fer)(x 2 - b 4 y 2 ) - 2pxcos 2 2fer 



2<ro sin 2 2br - —I 2 1 2 sin 26r 
yy 26 



Finally, the combination of equations (|6.19p , (|6.32|) and (|6.34|) leads to the 
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superposition formula 

1 / ,9/ sin6(ri — r 2 ) , 

y = ^O* - *i + b2 (yi + *)) ^['1%^ 

pi sin 26r 2 - p 2 sin 2bn ) 

U = Ml + W 2 + 2 — — ■ £ 

sm Zo(ti + t 2 

1 / 2ri sin 2 (2for 2 ) 2r 2 sin 2 (26r!) sin 2br\ sin 2&r 2 \ 2 (535) 

+ 62 Uin 2 (26(r 1 +r 2 )) + sin 2 (26(r 1 + r 2 )) ~ 6sin2&(r 1 + r 2 ) ) * ' 

T\ sin 26r 2 — t 2 sin 2br\ 

P = Pi+P2 + 2 , 2 . , , ; ; 4, 

(r sm2o(ri + r 2 ) 

2 2 T"j sin 26t 2 + r 2 sin 2&Ti 

q= -b Pi+b p 2 - 2 — ■ £, 

sm2o(ri + t 2 ) 

s = -n+r 2! r = ri+T 2 , where £ = i(x 2 - xi - fe 2 (yi + y 2 )). 

A A1 AO 

The restriction of V to the manifold Mi = {I = I = 0} gives 

W = { dui — pi dx\ + b 2 dyi, dpi — l/& 3 (tan6ri — brijdxi + Tidy%, (6.36) 
n dxi - b(bn + cot(6n)) dyi - b 2 dp x }. (6.37) 

This is a rank 3 Pfaffian system on a 5 manifold whose derived flag has dimen- 
sions [3, 2, 0]. The equivalence problem for such systems was analyzed in detail 
by Cartan [7], where it is established that the fundamental invariant for such 
systems is a certain rank 4 symmetric tensor T in two variables. For b = 0, this 
tensor vanishes while for b =^ we find that T is the 4-th symmetric power of 
a 1-form. In accordance with Cartan's result the symmetry algebra of W when 
b = is the 14 dimensional exceptional Lie algebra <? 2 and, indeed, it is not 
difficult to transform W to the canonical Pfaffian system for the Hilbert-Cartan 
equation z = (y ) . For b ^ the symmetry algebra of W is the 7 dimensional 
solvable Lie algebra with infinitesimal generators 

{d xl ,d yi ,d ux ,x 1 d Xl +yid yi + 2^8^ +p 1 d pi ,(x 1 - 6 2 yi)9 ui ~ d pi ,Y 1 ,Y 2 }, 

(6.38) 

where F 2 = and 

Y\ ={x\ + b 2 yi){bcot br\d Xl — -tanfrnc^ + 2pibcscbrid Ul + 2t\ csc26ti<9 Pi ) 



+ 2xiy%d ui + (yi - -^;)d Pl - d Tl . (6.39) 
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The PfafRan (|6.37p with b ^ may be transformed into Cartan [7], page 
170, equation (5'). 

Example 6.3. For our next example, let G be an n-parameter matrix group 
and, for the mapping (x,y) — > U(x,y) € G, consider the system of differential 
equations 

U X y = UJJ^Uy. (6.40) 

The general solution to these equations is well-known to be U (x, y) = A(x)B(y), 
with A(x), B(y) e G. In the case when U is a 1 x 1 matrix, this system reduces 
to the wave equation v xy = under the change of variable u = exp(u). We show 
how our integration method leads directly to the general solution and, in the 
process, we calculate the Vessiot algebra of (|6.40[) to be the Lie algebra of G. 
The Pfaffian system for (|OD|) is I = {9, 6\ 9 2 }, where 

Q = dU — U x dx — U y dy, 
9 1 = dU x - U xx dx - UJJ^Uy dy, (6.41) 

6 2 = dUy - UJJ^Uy dX - Uyy 6^. 

The first integrals for the singular systems are 

P=y, P=U- 1 U y , P=D y (P) = U- 1 U yv -U- 1 U y U- 1 U y , 
P = x, P= U x U-\ P = D X (P) = U^U- 1 ~ UJJ^UJJ- 1 , 

and our 0-adapted coframe for I is {0, dl ,dl , fj, dl ,dl , fj}, where 

fj = df - PdP = U^e 2 - If- 1 6 P, and 

rj = dP - PdP = e 1 ^- 1 - P e u- 1 . 

The structure equations are 

d® = dl 1 A(U fj + Q P) + dP A(fjU + PQ). (6.43) 

This coframe satisfies (14.13[) and is therefore 2-adapted. 

The next step is to eliminate either the dl A (I 9) or the dl A (0 / ) terms 
in (|6.43[) . By inspection, we see that the forms Qx = U^ 1 ® and Oy = QU^ 1 
provide us with the required 4-adapted coframes, with structure equations 

dO x = dP A fj + dP A {U^fj U) - O x A ® x + dP A {Q X P - P®x), 

dO Y = dP A (U fj U- 1 ) + dP Af) + <dy A Oy — dP A {QyP - P®y). 

(6.44) 
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Since the forms &x and 6y are Lie algebra valued, these structure equations 
show that the Vessiot algebra for (|6.40[) is the Lie algebra of G (Theorem 14.61) . 
The final coframe adaptation in Section 4.4 is given by 

6 = <d x + Pdl 1 = U^dU - XJ- l U x dx and 

(6.45) 

e = Y + I dr = dU u- 1 - U y U~ l dy, 
with structure equations 

dQ = dl 1 A {U- 1 ?) U)-QAQ and dO = dl 1 A (U fj U' 1 ) + Q A 6. (6.46) 

The form (|5 . 30[) are then found to be precisely the left and right invariant forms 
on G, that is, 

u} = U~ x dU and cj = dJ7C/" 1 , (6.47) 

so that the map p constructed in Theorem 15. H is simply p{x 1 y, U, . . .) = U. 

With respect to coordinates x,U\,U\ x ,U\ xx on the M\ = {I a = 0} and 
coordinates y, U% v , U2 Vy on the level set M2 = { 1 = 0}, the Pfaffian systems 

A v 

W and W are 

W = {dUi-Ui x dx, dUi x -U lxx dx} and W = {dU 2 -U 2 ydy 1 dU 2 y-U 2 yydy} 
and the superposition formula is 

U = U1U2, U X = Ul X U 2l Uy — U\U 2 y, U XX — Ul XX U 2l Uyy = U \U 2 yy . 

Example 6.4. It is an open problem to determine which scalar Darboux inte- 
grable equations admit generalizations wherein the dependent variable U takes 
values in an arbitrary non-commutative, finite dimensional algebra A. Here are 
two such examples which provide us with many Darboux integrable systems 
amenable to the methods presented in this paper. 

I- U Xy =(U X +l)U- 1 Uy 

II. U xy = U x {U-y)- 1 U y + U y (U-x)- 1 U x . 

The first integrals for the singular systems (excluding x and y) and general 
solutions are|f| 



3 We remark that the solution to II given by Vessiot 28 (equations C|» and (70), pages 6 
and 45) in the scalar case is incorrect. 
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I. f = (U x + I)U-\ I 2 = D x (I l ), T = U- x U yy , 
U = (F')- 1 {-F + G), 

II. I 1 = (U - x)- 1 [U xx U x l (U - x) - 2U X + I}, 
I 1 = (U-y)- 1 {U yy U; 1 (U-y)~2U y +I} : 
U = {xF 1 + V G'-F-G) (F' + G')- 1 , 

where F = F(x) and G = G(y) take values in A. For both systems I and II 
the Vessiot algebra is the tensor product of A with the Vessiot algebra for the 
corresponding scalar equation. We conjecture that all equations of Moutard type 
( |15j Volume II, page 250, equation 19) admit non-commutative generalizations. 

Example 6.5. Some of the simplest examples of Darboux integrable systems 
can be constructed by the coupling of a Darboux integrable scalar equation to 
a linear or Moutard- type equation. As examples, we give 



I. 


1l X y 


= e 2u , 


V X y 


= n(n -| 


- l)e 2u v, 




II. 


U X y 


6 Uy , 


^xy 


+ ((n- 


a)e u + au x )v x = 0, 




III. 


U X y 


6 Uy , 


^xy 


" e u v y - 


f (n+ l)u y v x + (n + 1)! 


e u u y = 0, 


IV. 


l^xy 


= e u u x 


^xy 


+ (el* 


- (nBe- v ) y + {n-2)B 


= 0, 



where B — e w u x and n is a positive integer. The system I appears in [19) (page 
116); systems II-IV do not seem to have appeared in the literature. 

For each of these systems the restricted Pfafhan systems W and W are jet 
spaces for two functions of a single variable (x or y). The Vessiot algebra for I is 
the semi-direct product of s[(2) and an Abelian Lie algebra of dimension 2n + 1, 
as determined by the (unique) (2n + l)-dimensional irreducible representation 
of s[(2). The infinitesimal action of the Vessiot group for I, restricted to W or 
W, is the action listed in [14] as number 27 (where now the variables x, y in 
[14] serve as the dependent variables for the jet spaces W and W). 

For II the Vessiot algebra is a semi-direct product of the 2-dimensional 
solvable algebra with an (n+ l)-dimcnsional Abelian algebra. The infinitesimal 
action of the Vessiot group, restricted to W or W, is number 24 in [14]. The 
infinitesimal Vessiot groups for III and IV have dimensions n + 3 and n + 4 and 
coincide, respectively, with numbers 25 and 26 in [14] , 
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For n = 1 the general solutions to these systems are 

_1 FjG[ _ n F 2 - G 2 F g G' 2 

11 M = h ^ ^-^(Fa-Ga-^-Gx)^; 

111 u = ln G7^' w =(f^-?X^M- lll(Gl - Fl) ' 

IV. u = In 





F X 






Gi- 


Fx 







Fx-Gx' 

= 1 f (Gj^-G^-G^Fx-Gx))^ 

Vi (F 2 - G 2 ) - F^ (F - d)) (F 1 - G x )>- 

The general solutions for arbitrary n can be obtained in closed compact form 
by the method of Laplace. 

In addition, any non-linear Darboux integrable system can be coupled to its 
formal linearization to obtain another Darboux integrable system. For example, 
if we prolong the partial differential equations 

V. 3u xx u yy + 1 = 0, v xx —v yy = 

u yy 

to order 3 in the derivatives of u, we obtain a rank 8 Pfaffian system on a 14- 
dimcnsional manifold which is Darboux integrable, with 4 first integrals for each 
associated singular Pfaffian system. 

The Vessiot algebra is Abelian and of dimension 6. The two Lie algebras of 
vector fields dual to the forms 6 and 9 for the 5-adapted coframe coincide and 
are given by 

{<9y, Out ~t~ &Uj. ) i ^ 0y x : UyOy -\-U X yOy x -\-UyyOy y -\-U X yyOv j . y -\~Uyyyd-u yy ^ . 

In accordance with Remark 14.101 these vector fields are also infinitesimal sym- 
metries for V. In terms of the arbitrary functions <p(a) and ip((3) appearing in 
the general solution to 3u xx u yy + 1 = (Goursat([U Vol. 2, page 130), the 
general solution for v is given implicitly as 

x = l^ZJ-> y=\{p- a ){<p" +r)+4>' *=^> (6-48) 

, — F \ C ^ ~ r + 03 " a)<t> '" F> f- ~ t. + {(3 - a)r ' G> (6 49) 
where F = F(a) and G = G{0). 
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Finally, we remark that the Laplace transformation (not the integral one, 
Darboux [9] , Forsyth [13], pages 45-59.) can be applied to the linear compo- 
nents of any of the above systems to obtain new Darboux integrability systems. 
As well, the linear component in any of these systems can be replaced by their 
formal adjoint to obtain yet other Darboux integrable systems. 

Example 6.6. Although we are unaware of an explicit general proof it is gen- 
erally acknowledged that the Toda lattice systems (see, for example, [T5], [22] ) 
are Darboux integrable. In this example we shall check that the B2 Toda lattice 
equations 

u xy = 2e u -2e v , v xy = -e u + 2e v (6.50) 

are Darboux integrable and find the closed-form, general solution. We use this 
example to illustrate a slightly different computational approach, one based 
upon the symmetry reduction interpretation of the superposition formula given 
in Sections 3 and 5.3. 

The canonical Pfaffian system for (|6.50|) satisfies our definition of Darboux 
integrable upon prolongation to 4-th order, that is, as a rank 14 Pfaffian system 
/ on a 20 dimensional manifold. The diffcomorphism x <-> y interchanges the 
singular Pfaffian systems V and V. The first integrals for the singular Pfaffian 

A 

system V (containing dx) are 

01 02 2 1 1 , 1, « , -. 02 £4 „ C3 n 

I = x, I =v xx + -u xx - -u x v x - -u x - -v x , I = D X I , I = D X I , and 
3 3 o 3 



I ^xxxx ~t~ 2v xxxx 2v x V xxx %L XX 1L X 2%L X 1L XX V X ~t~ U x 4" ^xx~^~ 



1 A 1 

2' 



"^222 ^"2 ^"3 

~^V X U X V x U X x 'q^^xx'^'x ^xx^xx H~ "T^x^x ^x^xxx- 

Let W be the restriction of V (or, equivalently, I) to Mi — {/ = 0}. We 

A 

find that W is a rank 12 Pfaffian system on a 15 manifold. The derived flag 
of W has dimensions [12, 10, 8, 6, 4, 2, 1, 0] while the dimensions of the space of 
Cauchy characteristics for these derived Pfaffian sytems are [0, 2, 4, 6, 8, 10, 12, 
15]. By using the invariants of these Cauchy characteristics as new coordinates, 
we are able to write W in the canonical form 

W ~ { dui — ii2dui, dii2 — U2dui, dui — u^dx, du' 1 ~du'{dx 1 ...}. (6.51) 

Here there are 7 contact forms for U2 and 5 for u\. In these coordinates the 
integral manifolds of W are given by ui = Fi(x), — F[{x), . . .and U2 = 
F 2 (F 1 {x)), ii 2 = (F 2 )(Fi(z)), .... 
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The 10 dimensional infinitesimal Vessiot group, restricted to Mi, now takes a 
remarkably simple and well-known form - it is the infinitesimal conformal group 
o(3,2) acting on the 3-dimensional space (ui, 112,112) by contact transforms. 
(See, for example, Olver [21] page 473.) Explicitly, the generating functions for 
the infinitesimal action of the Vessiot group on M\ are 

Q = [U2; — «2 + Ul"2, -Uiul + 2u 2 U2; 2 ' ""2' 2^2' U l> 

_ 2 U i"2 _ 2u 2 + 2uiu 2 U2, ~2u\U2 + u\u 2 , 2 U \\- (6.52) 

A 

To obtain the infinitesimal generator X q corresponding to a function q £ Q, first 
construct the vector field X° = -qu 2 d Ul +(q— u 2 qu 2 )d U2 + (q u + U2<lu 2 )du 2 and 
then prolong X q to the vector field X q on M\ by requiring it to be a symmetry 
of W. We remark that the basis for o(3, 2) so obtained is the canonical Chevalley 
basis in the sense that the first 2 vectors define the Cartan subalgebra, the next 
4 correspond to the positive roots, and the last 4 to the negative roots. 

By Theorem 15.131 the superposition formula for the B 2 Toda lattice can 
therefore be constructed from the joint invariants for the diagonal action of the 
conformal algebra (3, 2) on M1XM2. We use coordinates [y, Vi, v{, v'[ . . . , V2, i>2, #2> • 
on M.2- To compactly describe these joint invariants we first calculate the joint 
differential invariants in the variables 

{ux, u[, vi, v[, ii 2 , 1)2, u 2 , v 2 , U2, V2} 

for the 7 dimensional subalgebra of o(3,2) generated by Qi, Q2, Qa, Q5, Qe, 
Q7, Qiq. These are 

Jl = v[(u 2 ) 1/3 (v2) 2/3 /(u2 ~ V 2 ), J2 = u[(u 2 ) 2 ^(v 2 ) 1/3 /(u2 ~ V 2 ), 

J 3 = (ti 2 «i - il 2 ui +u 2 ~ V2)(u 2 ) 1/3 (v 2 ) 2/3 /{v2 - u 2 ) 2 , 

J 4 = (v 2 ui - v 2 Vi +v 2 - U2)(u 2 ) 2/3 (v 2 ) 1/3 /(v 2 - u 2 ) 2 , and (6.53) 

•h = —(U2V2V 2 — 2U2V2V\U\ + U2V2u\ — 2u 2 «2Ul + 2lL2U2V\ + 2it 2 M2 

- 2u 2 V2 - 2v 2 V2V\ + 2v 2 V2Ui + 2v 2 V2 - 2v 2 u 2 (6.54) 
+ ill - 2u 2 V2 + vl)u2V2/(2{v2 - u 2 ) A ). 

Then, in terms of these partial invariants the (lowest) order joint differential 
invariants for o(3, 2) are 

JiJ 2 (J 3 J 4 -2J 5 ) 2 JiJ 2 (J 3 J4- J 5 ) 

Kl = — m — — and K2 = — m — ^n~u~ 6 - 55 

(J3J4-J5) (J 3 J 4 -2J 5 )^ 
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and the solutions to the B 2 Toda lattice are 

u = \n(K 1 /4) and v = \n(2K 2 ), where (6.56) 
ui = F^x), u 2 = F 2 {F 1 {x)), Vl = Gi(tf), v 2 = G 2 {G 1 {y)). (6.57) 

It is hoped that a more transparent representation of these solutions, similar to 
that available for the A n Toda lattice will be be obtained. 

Example 6.7. Let P denote the 2-dimensional Minkowski plane with metric 
dx Q dy and let N be a pseudo-Riemannian manifold with metric g. A map- 
ping ip : P — > iV which is a solution to the Euler-Lagrange equations for the 
Lagrangian 

L = 9[ %% )dx ^ dy (6 - 58) 
is called a wave map. There are precisely two inequivalent, non-flat metrics (up 
to constant scaling) in 2 dimensions, namely 

.91 = 1 ^ (du 2 + dv 2 ) and g 2 = -^—^(du 2 + dv 2 ) (6.59) 

which define Darboux integrable wave maps at the 2-jet level (that is, without 
prolongation). Surprizingly, these metrics are not constant curvature. It is not 
difficult to check that under the change of coordinates 



x = x — y, t = x + y, 9 = arctan(V e u — 1), \ — v /% (6.60) 

the differential equations (4.13) in [3] become the wave map equations for the 
metric g 2 . The Vessiot algebras for the wave map equations for g\ and 172 are 
st(2) x R and so(3) x R respectively. 
The wave map equations for g\ are 

V X Vy ~ U x Uy U X Vy + UyV x . . 

u xv = -, v xv = — . (6.61) 

v 2e u + 2 ' y 2e u + 2 y ' 

The standard encoding of these equations as a Pfaffian system results in a rank 
6 Pfaffian system on a 12 manifold. There are 4 first integrals for each singular 
Pfaffian system - for the singular Pfaffian system V containing dx the first 

1 n f2 eU ( u l + v l) ?3 „ f2 , 

integrals are I — x, I = — — — , I = D X I , and 

f + e u 

f4 v x u x (1 + 2e u )v x 

1 = ..2 I „.2 U ^ ~ ..2 I Z5 V ** 2 + 2e u ■ ( 6 - 62 ) 



M X 1 X 
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After considerable computation, the superposition is obtained and we find the 
general solution, in terms of the four arbitrary functions F\(x), F 2 (x), Gi(y), 
G 2 {y) to be 

2e u = -1 + v 7 ! + AVl + B 2 + AB sin(A) (6.63) 

and 

v = F 1 {x)-F 2 {x) + G 1 {y)-G 2 {y) (6.64) 



arctan M^l±^ + arctan ^ + B 



+ arctan 

+ arctan 
where 



AB' cos(A) + G' 2 B 2 Vl + A 2 Vl + B 2 + G' 2 AB(1 + B 2 ) sin(A) ' 
G^ABVTTB 2 cos(A) - BB'Vl + A 2 - AB'VTTB 2 sin(A) 

' A'B cos(A) - F^A 2 Vl + A 2 Vl + B 2 - F^AB(1 + A 2 ) sin(A) \ 
F^ABVTTA 2 cos(A) + AA'Vl + B 2 + A'BVTTA 2 sin(A) i 



.l(.r)=,/('§V-2§ > B(») = W f - 2^, A = F 2 -G 2 . (6.65) 



^2/ ^2" VV^V G 2 



Note that 



F{(x) = F^x)(l + ./l + A(x) 2 ), G' 1 {y) = G' 2 {y){l + ^/TTBW). (6.66) 
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Example 6.8. We turn now to some simple examples of overdetermined sys- 
tems for a single unknown function of 3 independent variables, beginning with 
the system 

U XZ = UU X , Uy Z = UUy. (6.67) 

The structure equations for the canonical encoding of this system as a rank 4 
Pfaffian system I = {a 1 , a 2 , a 3 , a 4 } on an 11-manifold are (modulo I), da 1 = 0, 

da 2 = n 1 At} 3 +tt 2 Att 4 , da 3 = tt 1 A n 4 + it 2 A tt 5 , da 4 = it 1 A tt 2 , (6.68) 

where tt 1 = dx, tt 2 — dy, tt 1 = dz, 

tt 3 = du xx - (u xx u + u 2 x )dz, tt 4 = du xy - (u xy u + u y u x ) dz, 

(6.69) 

tt° = duyy — (uyyU + u y ) dz , tt 2 = du g z — (u z + u 2 ) (u x dx + Uy dy) . 

The first integrals for the singular Pfaffian systems V = I U {tt 1 , . . . , 7T 5 } and 
V = I U {tt 1 , tt 2 } are I 1 — x, I 2 — y, I 1 = z, 

i = — , 7 = D x l , 1 — Dyl , 1 —u z — u 2, 1 = D z l . (6.70) 

u x 

The form tt 3 is not in V°° + V and therefore (|6.67p is not Darboux integrable 
on the 2-jet. The prolongation of (|6.67|) defines a decomposable rank 8 Pfaffian 
system /W = {a 1 , . . . ,a 8 } on a 16 dimensional manifold. In addition to the 
first integrals (|6.70[) . we now also have 

f = D 2 J 3 , V = D xy l\ ? = Dyl 3 , f = ^ - , 7 4 = DIP (6.71) 

^a; ^U x 

and the conditions (| 1 ,8p (| 1 . 10[) for Darboux integrability of are now satisfied. 
A 1-adapted coframe is a 1 = dl 1 , a 2 = dl , a 6 — dP , a 7 — dl , a 8 = dI S , 



f7 J = 


= df, 


a 1 = dl\ a 2 = dT, 


~1 
'/ 


i 

u x 


3 U V 2 >f3 r 4 ,fl f 5 ^f2 

a §- a = ai — 1 dl — 1 dl , 

u x 


^2 

V 


1 

Ux 


< u% 




1 

U x 


^ Uy^ U^y^ -2UyU X y + Uyy U Xa2=dF f7 ^ f ^ 

u x U x U x 


v 1 


= a 4 


1 ,r2 f3 ,yi -2 8 4 1 .f3 f4 ,ri 
— ua — dl — 1 dl , i] —a—ua — u z a = dl — 1 dl , 


e 1 


= a 1 


— du — u x dx — u y dy — u z dz, 


o 2 


2 

= a 


= du x — u xx dx — u xy dy — uu x dz, 


e 3 


= a 3 


du xx u xxx dx u XX y dy [uu xx -t- u x ) dz. 
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This frame is in fact 2-adapted. We next calculate U = U U {Xi, X2, X3} 
and U {Qc) = UU {Yi,Y 2 , Y$}, where 

v ?4 a . ^4 Ag 3< x i + 4u xx u x I -2u x I I 
Xi = -u x I d e i - (u xx I + u x I )d 9 2 — - d„ 



lUr. 



X 2 = -u x dgs, X 3 = -u x d e i - u xx d g 



3u 2 xx + 2u 2 x f 
2u T 



d e3 , Y 1 = -d e i, (6.72) 



Y 2 = udgi + u x d 2 + u xx d 8 3, Y 3 = ( — — + I 2 )d e i - uu x d S 2 - (u 2 x + uu xx )d 6 3. 

For the computations of Section 4.3 we use the base point defined by setting 
u = 0, u x = 1, u xx = 0, I = 1, and all other first integrals (|6.70p - (|6.71j) to 0, 
The matrices (|4.33j) and (|4.36|) are then given by 





' 1 


2/ 4 ? 


f 


p = 


F 


J 6 


F 





1 










-J 9 


1 



and i? 



0-10 
J 2 1 
10 



(6.73) 



From these matrices we calculate the 4 adapted coframes 







1 




u x 


— ( 




u x 




V 2 


2 




"SI* 


+ - 















l xx a 2 



9 2 



U T + uu 



xx a 2 



W) 



l xx a 2 



u(2u 2 + uu a 
2^! 



(6.74) 



2ul 



The Vessiot algebra is sl(2). Because this algebra is semi-simple, we can use 
Case I of Section 4.4, Theorem 15.11 and (|6.74p to directly determine the map 
p : M -> Aut(st(2) as 



p(u,u x ,u xx ) = A 



r 2 

U x UUxx 


u 


2 

u x 


u x 




1 


u l 


u x 


u(uu xx - 2u 2 x ) 


u 2 


[ 2u 2 x 


2u x 



u xx {uu xx - 2u 2 x )~ 



2 f x 

u xx 

(uu xx - 2u x ) 2 
4ul 



To obtain the superposition formula we introduce local coordinates (z, w, 

nil, w xx , w z , w zz , w zzz ) for Mi and ( 

•^1 Hi V) Vx-, V Zl VxXl ^Xy-, Vyy-, Vyzi ^XXXi 



Examples 



77 



v X xy, v xyy , Vyyy, V yyz ) for M% . The inclusions L\ : M\ — > M and 62 '■ M2 — > M 
are fixed by 

li(wi) = ui, ti(/ a )=0, i2{vi) = ui 1 L 2 '{I a ) = 0. 
The superposition formula is then found by solving the equations 

i* 1 {I a )=I a 1 £,{i a )=I a , p(u,u x ,u xx ) = p{w, w Xl w xx ) -p{v,v x ,v xx ) (6.76) 
for the coordinates of M . We find that 

u = w - . (6.77 

—2w x + vw xx 

Finally, we calculate the integral manifolds for W and W. It is immediate 
that W is the canonical Pfaffian system on J 3 (R, R 2 ) and hence the integral 
manifolds are defined by v = V{x, y). The last non-zero form in the derived flag 
for 

/\ 

W = { dw — w z dz, dw z — w zz dz, dw zz — w zzz dz, dw x — w x wdz, 
dw xx + (-wl - ww xx )dz } 

yields the Pfaffian equation 



dw xx — ^^-dw x — w 1 dz = 



from which it follows that w x — G'(z) and w xx = G(z)G' (z). One of the 
remaining equations in W then gives w = G"(z)/G'(z). On replacing V(x,y) 
by —2/F(x,y), the superposition formula (|6.77p becomes 

G"(z) 2G'(z) 
u{x ^ Z) =mz)- F{ X ,y) + G{zY (6J8) 

which gives the general solution to (|6.67p . 

We continue this example by considering two variations on (|6.67p . First we 
observe that to (|6.67p we may add any equation of the form 

p( x y — UyUxx ~ UxUx v u v u xy ~~ u x u tyy . _ Q ,g 

U x 1l x 

to obtain a rank 4 Pfaffian system I — {a 1 , a 2 , a 3 , a 4 } on a 10 dimensional 
manifold. The structure equations become ( mod /) da 1 = 0, 

da 2 = n 1 An 2 , da 3 = tt 3 A # 4 , da 4 = n 1 A tt 2 (6.80) 
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and hence all such systems are involutive with Cartan character s\ = 3. 
For example, consider the system of 3 equations 

^yHxy ^x^yy Oj ^xz ^^x: ^yz Utly (6.81) 

The foregoing calculations can be repeated, almost without modification, to 
arrive at the same superposition formula (|6.77|) - the only difference is that now 
W is the (prolonged) canonical Pfaffian system for the equation v y v xy — v x v vv = 
0, an equation which is itself Darboux integrable. Thus, in more complicated 
situations, the method of Darboux, can be used to integrate the Pfaffian systems 
W and W and the superposition formula for the original system is given by a 
composition of superposition formulas. In the case of the present example, 
the calculation of the first integrals for W reveals that this system is contact 
equivalent to the wave equation (via X = x, Y = v, V = y, Vx = —v x /v v , 
Vy = l/^i/) and leads to the parametric solution 

G"(z) 2G'{z) 
x = a, y = f(a)+g(r), u = — (6.82) 

Our second variation of (|6.67p is obtained by the differential substitution 
u x = exp(w). This leads to the equations 

v xz = exp(v) and v yzz = v yz v z . (6.83) 

It is surprising that the canonical Pfaffian system for these equations (obtained 
by the restriction of the contact ideal on J 3 (R, R 2 )) does not define a decom- 
posable Pfaffian system. The following theorem resolves this difficulty. 

Theorem 6.9. The system of differential equations 

u xz = F(x, y, z, u, u z ,u yz ,u zz ), 

^■IJZZ C(*£? ?/, Z, Zi, U X7 tiy, U Z: Uy Z: U ZZ7 Uyy Z7 11 Z Z Z ) 

determines the rank 4 Pfaffian system 

a 1 = du — u x dx — u y dy — u z dz, a 2 = du z — Fdx — u yz dy — u zz dz, 

a 3 = du yz — D y (F)dx — u yyz dy — Gdz 7 (6.85) 

a A = du zz — D z (F)dx — Gdy — u zzz dz, 

on an 11 dimensional manifold. If the compatiblity conditions for (|6.84p hold, 
then this Pfaffian system is decomposable and involutive with Cartan characters 
s\ — 1 and S2 = 1. 



(6.84) 
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We use Theorem 16.91 to write (|6.83l) as a rank 4 Pfaffian system on an 1 1- 
manifold. The prolongation of this system is Darboux integrable and calcula- 
tions, virtually identical to those provided for (|6.67|) . lead directly to the general 
solution 
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